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ABSTRACT 

We present the third Fermi Large Area Telescope source catalog (3FGL) of 
sources in the 100 MeV-300 GeV range. Based on the first four years of science 
data from the Fermi Gamma-ray Space Telescope mission, it is the deepest yet in 
this energy range. Relative to the 2FGL catalog, the 3FGL catalog incorporates 
twice as much data as well as a number of analysis improvements, including 
improved calibrations at the event reconstruction level, an updated model for 
Galactic diffuse 7 -ray emission, a rehned procedure for source detection, and im¬ 
proved methods for associating LAT sources with potential counterparts at other 
wavelengths. The 3FGL catalog includes 3033 sources above da signihcance, 
with source location regions, spectral properties, and monthly light curves for 
each. Of these, 78 are flagged as potentially being due to imperfections in the 
model for Galactic diffuse emission. Twenty-five sources are modeled explicitly 
as spatially extended, and overall 238 sources are considered as identified based 
on angular extent or correlated variability (periodic or otherwise) observed at 
other wavelengths. For 1010 sources we have not found plausible counterparts at 
other wavelengths. More than 1100 of the identihed or associated sources are ac¬ 
tive galaxies of the blazar class; several other classes of non-blazar active galaxies 
are also represented in the 3FGL. Pulsars represent the largest Galactic source 
class. From source counts of Galactic sources we estimate the contribution of 
unresolved sources to the Galactic diffuse emission is ~3% at 1 GeV. 

Subject headings: Gamma rays: general — surveys — catalogs 


1. Introduction 

This paper presents a catalog of high-energy 7 -ray sources detected in the first four 
years of the Fermi Gamma-ray Space Telescope mission by the Large Area Telescope (LAT). 
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It is the successor to the LAT Bright Source List (hereafter OFGL. lAbdo et ahll2009d[) . the 
First Fermz LAT (IFGL, Abdo et ah 2010d) catalog, and the Second Fermi LAT (2FGL, 
Nolan et ah 2012h catalog, which were based on 3 months, 11 months, and 2 years of flight 
data, respectively. The 3FGL catalog both succeeds and complement s the First Fermi LAT 
Gatalog of Sources Above 10 GeV (IFHL, lAckermann et aLll2013al) . which was based on 
3 years of flight data but considered only sources detected above 10 GeV. The new 3FGL 
catalog is the deepest yet in the 100 MeV-300 GeV energy range. The result of a dedicated 
effort for studying the Active Galactic Nucl ei (AGN) population i n the 3FGL catalog is 
published in an accompanying paper fOLAG. lAckermann et al.ll2015l) . 


We have implemented a number of analysis rehnements for the 3FGL catalog: 


1. Pass 7 reprocessed datal^ are now used (§ 12.21) . The principal difference relative to the 
original Pass 7 data used for 2FGL is improved angular resolution above 3 GeV. In 
addition, systematics of the instrument response functions (IRFs) are better charac¬ 
terized, and smaller. 


2. This catalog employs a new model of the diffuse Galactic and isotropic emissions, 
developed for the 3FGL analysis (§ 12.31) . The model has improved hdelity to the 
observations, especially for regions where the diffuse emission cannot be described using 
a spatial template derived from observations at other wavelengths. In addition, the 
accuracy of the model is improved toward bright star-forming regions and at energies 
above 40 G eV generally. The development of this model is d escribed in a separate 
publication ffGasandiian fc the Fermi LAT Gollaboration 2015 ). 


3. We explicitly model twenty-hve sources as extended emission regions (Ij 13.4|) . up from 
twelve in 2FGL. Each has an angular extent measured with LAT data. Taking into 
account the hnite sizes of the sources allows for more accurate flux and spectrum 
measurements for the extended sources as well as for nearby point sources. 


4. We have further rehned the method for characterizing and localizing source ‘seeds’ 
evaluated for inclusion in the catalog (§ 13.ip . The improvements in this regard are 
most marked at low Galactic latitudes, where an iterative approach to Ending seeds 
has improved the sensitivity of the catalog in the Galactic plane. 


5. For studying the associations of LAT sources with counterparts at other wavelengths, 
we have updated several of the catalogs used for counterpart searches, and correspond¬ 
ingly recalibrated the association procedure. 


^See http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Pass7REP_usage.html 
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The exposure of the LAT is fairly uniform across the sky, but the brightness of the diffuse 
backgrounds, and hence the sensitivity for source detection, depends strongly on direction. 
As for previous LAT source catalogs, for the 3FGL catalog sources are included based on 
the statistical signihcance of their detection considered over the entire time period of the 
analysis. For this reason the 3FGL catalog is not a comprehensive catalog of transient 7 -ray 
sources, however the catalog does include light curves on a monthly time scale for sources 
that meet the criteria for inclusion. 

In § [2] we describe the LAT and the models for the diffuse backgrounds, celestial and 
otherwise. Section |3] describes how the catalog is constructed, with emphasis on what has 
changed since the analysis for the 2FGL catalog. The 3FGL catalog itself is presented 
in § m along with a comparison to previous LAT catalogs. We discuss associations and 
identideations in § O and Galactic source counts in § [HI The conclusions are presented in 
§ [71 We provide appendices with technical details of the analysis and of the format of the 
electronic version of the catalog. 


2. Instrument and Background 
2.1. The Large Area Telescope 

The LAT detects 7 rays in the energy range 20 MeV to more than 300 GeV, measuring 
their arrival times, energies, and directions. The LAT is also an efficient detector of the 
intense background of charged particles from cosmic rays and trapped radiation at the orbit 
of the Fermi satellite. Accounting for 7 rays lost in hltering charged particles from the 
data, the effective collecting area is ~6500 cm^ at 1 GeV (for the P7REP_SOURGE_V15 
event selection used here; see below). The live time is nearly 76%, limited primarily by 
interruptions of data taking when Fermi is passing through the South Atlantic Anomaly 
(~13%) and readout dead-time fraction (~9%). The held of view of the LAT is 2.4 sr at 
1 GeV. The per-photon angular resolution (point-spread function, PSF, 68 % containment) 
is ~5° at 100 MeV, decreasing to 0?8 at 1 GeV (averaged over the acceptance of the LAT), 
varying with energy approximately as 77“°'® and asymptoting at ~0?2 above 20 GeV. The 
tracking section of the LAT has 36 layers of silicon strip detectors interleaved with 16 layers 
of tungsten foil (12 thin layers, 0.03 radiation length, at the top or Front of the instrument, 
followed by 4 thick layers, 0.18 radiation length, in the Racfcsection). The silicon strips track 
charged particles, and the tungsten foils facilitate conversion of 7 rays to positron-electron 
pairs. Beneath the tracker is a calorimeter composed of an 8 -layer array of Gsl crystals 
(~8.5 total radiation lengths) to determine the 7 -ray energy. A segmented charged-particle 
anticoincidence detector (plastic scintillators read out by photomultiplier tubes) around the 


tracker is used to reject chargee 
L AT is provide c 


in 


in 


Atwood et al 


- partic le background events. More information about the 
fl2009h. a nd the in-fli ght calibration of the LAT is described 


Abdo et ahl f 2009d) . I Acker man n et a 1.1 f2012 d) . and lAckermann et ahl (2012a). 


2.2. The LAT Data 


The data for the 3FGL catalog were taken during the period 2008 August 4 (15:43 
UTC)-2012 July 31 (22:46 UTC) covering close to four years. They are the public data 
available from the Fermi Science Support Ce nter (FSSC). Intervals a round bright GRBs 
(080916C, 090510, 090902B, 090926A, 110731A: lArkermann et allboiObh were excised. Solar 
flares became relatively frequent in 2011-12 (close to solar maximum) and were excised as 
well whenever they were bright enough to be detected over a month at the 3a level. Solar 
flares last much longer than GRBs, so we were attentive not to reject too much time. Since 
the 7 -ray emission is localized on the Sun, we kept intervals during which the Sun was at 
least 3° below the Earth hmh@ even during solar flares. The solar flares were detected over 
3-hour intervals so the corresponding Good Time Intervals (GTI) are aligned to those 3- 
hour marks. Overall about two days were excised due to solar flares. In order to reduce 
the contamination from the y-ray-bright Earth limb (Ij I2.3.4|) . times when the rocking angle 
of the spacecraft was larger than 52°, and events with zenith angles larger than 100°, were 
excised as well. The precise time intervals corresponding to selected events are recorded in 
the GTI extension of the 3FGL catalog (FITS version, App. iBj) . 

The rocking angle remained set at 50° after September 2009 (it was 35° for about 80% 
of the time before thallfl). With the larger rocking angle the orbital plane is further off axis 
with the result that the survey is slightly non-uniform. The maximum exposure is reached 
at the North celestial pole. At 1 GeV it is 60% larger than the minimum exposure, which is 
reached at the celestial equator (Figure [T]). 

In parallel with accumulating new data, developments on the instrument analysis side 
fiBregeon et al.ll2013l) led to reprocessing all LAT data with new calibration constants, re¬ 
sulting in the Pass 7 reprocessed data that were used for 3FGL0. The main advantage for the 
source catalog is that the reprocessing improved the PSF above 10 GeV by ~30%, improving 


^This selection in FTOOLS notation is ANGSEP(RAJUN,DEC_SUN,RA_ZENITH,DEC_ZENITH) > 115. 

^See the LAT survey-mode history at http: //f ermi. gsf c .nasa.gov/ssc/observations/types/allsky/, 

"‘Details about the performance of the LAT are available at 
http://www.slac.Stanford.edu/exp/glast/groups/canda/lat_Performance.htm. 
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16 18 20 22 


Fig. 1.— The exposure at 1 GeV in Galactic coordinates and Hammer-Aitoff projection 
for the four-year period analyzed for the 3FGL catalog. The units are equivalent on-axis 
observing time (in Ms). 


the localization of hard sources (§ 13.ip . We used the Source class event selection. 

The lower bound of the energy range was left at 100 MeV, but the upper bound was 
raised to 300 GeV as in the IFHL catalog. Thi s is because as th e source-to-background 
ratio decreases, the sensitivity curve (Figure 18 of lAbdo et al.ll2010di IFGL) shifts to higher 
energies. 


2.3. Model for the Diffuse Gamma-Ray Background 

Models for the diffuse y-ray backgrounds were updated for the 3FGL analysis, taking 
into account the new IRFs for Pass 7 reprocessed data and the improved statistics available 
with a 4-year data set, and also applying some rehnements in the procedure for evaluating the 
models. The primary components of the diffuse backgrounds are the diffuse y-ray emission 
of the Milky Way and the approximately isotropic background consisting of emission from 
sub-threshold celestial sources plus residual charged particles misclassihed as y rays. In 
addition, we treat the ‘passive’ emission of the Sun and Moon from cosmic-ray interactions 
with the solar atmosphere, solar radiation held, and the lunar lithosphere, as effectively a 
diffuse component, because the Sun and Moon move across the sky. The residual Earth limb 
emission after the zenith angle selection (§ 12.2p is also treated as effectively diffuse. Each 






















component of the background model for the 3FGL analysis is described in more detail below. 


2.3.1. Diffuse emission of the Milky Way 


The diffuse 7 -ray emission of the Milky Way originates in cosmic-ray interactions with 
interstellar gas and radiation. As for 2FGL, for any given energy the model is primarily a 
linear combination of template maps derived from GO and H l line survey data plus infrared 
maps of interstellar dust, which trace interstellar gas and in the model represent the 7 - 
ray emission from pion decay and bremsstrahlung. In addition we include in the model a 
template representing the intensity of emission from inverse Gompton scattering of cosmic- 
ray electrons on t he interstellar radiation field . This component was calculated using the 
GALPROP cod^ ( Moskalenko fc Strong 1998). 


For the 3FGL analysis we have made several improvements relativ e to 2FGL in modeling _ 

the di ffuse emission. The development of this new model is described bv iGasandiian fc the Fermi LAT Golla 
( 2015 1^. Here we briefly summarize the improvements. For 3FGL the representation of the 
gas traced uniquely by the infrared maps was improved in the vicinity of massive star-forming 
regions. The overall model was £t to the LAT data iteratively, taking into account a prelim¬ 
inary version of the 3FGL source list, for the energy range 50 MeV-50 GeV. Relative to the 
2FGL model for Galactic diffuse emission, we have improved the extrapolation to lower and 
higher energies using the energy dependence of the 7 -ray emissivity function. We developed 
a new procedure to account for the structured celestial 7 -ray emission that could not be fit 
using templates derived from observations at other wavelengths. This residual component, 
also fit iteratively, was derived by deconvolving the residuals to take into account the effects 
of the PSF, filtering the result to reduce statistical fluctuations (removing structures on an¬ 
gular scales smaller than ~2°). The spectrum was modeled as inverse Gompton emission 
from a population of cosmic-ray electrons with a spectral break at 965 MeV. 


2.3.2. Isotropic background 

The isotropic diffuse background was derived from all-sky fits of the four-year data set 
using the Galactic diffuse emission model described above and a preliminary version of the 
3FGL source list. The diffuse background includes charged particles misclassihed as 7 rays. 


®See http://galprop.stanf ord.edu. 


®The model is available as gll_iem_v05_revl. f it from the FSSC. 
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We implicitly assume that the acceptance for these residual charged particles is the same as 
for 7 rays in treating these diffuse background components together. For the analysis we 
derived the contributions to the isotropic background separately for Front-converting and 
Facfc-coaverting events. They are available as iso_source_xxx_v05.txt from the FSSC, 
where xxx is front or back. 


2.3.3. Solar and lunar template 


The quiescent Sun and the Moon are fairly bright y-ray sources flAbdo et ahl l2011b 


2012b|). The Sun moves in the ecliptic but the solar y-ray emission is extended because 
of cosmic-ray interactions with the solar radiation held; detectable emission from inverse 
Compton scattering of cosmic-ray e l ectron s on the radiation held of the Sun extends several 
degrees from the Sun flAbdo et al.N2011bl) . The Moon is not an extended source in this 
way but the lunar orbit is inclined somewhat relative to the ecliptic and the Moon moves 
through a larger fraction of the sky than the Sun. Averaged over time, the y-ray emission 
from the Sun and Moon trace a region around the ecliptic. We used models of their observed 
emission together with calculations of their motions and of the exposure of the observations 
by the LAT to ma ke templates for the equ ivalent dihuse component for the 3FGL analysis 
using gtsuntemp fjjohannesson et ah 2013 1. For the light curves (§ 13.Op we evaluated the 
equivalent diffuse components for the corresponding time intervals. 


2 . 3 . 4 . Residual Earth limb template 


The limb of the Earth is an intense source of y rays from cosmic-ray collisions with 
the upper atmosphere flAbdo et al.ll2009all . At the ~565 km altitude of the (nearly-circular) 
orbit of the LAT, the limb is ~112° from the zenith. During survey-mode observations, 
which predominated in the hrst four years of the Fermi mission, the spacecraft was rocked 
toward the northern and southern orbital poles on alternate ~90 minute orbits. With these 
attitudes, the edge of the LAT held of view closest to the orbital poles generally subtended 
part of the Earth limb. As described in § 12.21 we limited the data selection and exposure 
calculations to zenith angles less than 1 00°. Because the Earth l imb emission is so intense, 
and the tails of the LAT PSF are long flAckermann et al.ll2012al) . a residual component of 
limb emission remained in the data. Over the course of a precession period of the orbit 
(~53 d), the residual glow hlls out large ‘caps’ around the celestial poles, with the angular 
radius determined by the sum of the orbital incl ination (25?6) and the angular distance of the 
zenith angle limit from the orbital pole (10°). ICasandiian fc the Fermi LAT Collaborati^ 
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fl2015h describe how the map and spectrum of the residual component were derived. The 
spectrum is well modeled as a steep power law in energy with index 4.25. This is steep 
enough that the residual Earth limb emission contributes signihcantly only below 300 MeV. 


3. Construction of the Catalog 

The procedure used to construct the 3FGL catalog has a number of improvements 
relative to what was implemented for the 2FGL catalog. In this section we review the 
procedure, with an emphasis on what is being done differently. The signihcances (§ 13.2p . 
spectral parameters (§ 13.31) and fluxes (§ 13.5p of all catalog sources were obtained using 
the standard pyLikelihood framework (Python analog of gtlike) in the LAT Science Tool^ 
(version v9r32p5). The localization procedure f G3.1l) . which relies on pointlike, provided the 
source positions, the starting point for the spectral htting, and a comparison for estimating 
the reliability of the results (§ 13.7.41) . Throughout the text we use the Test Statistic TS = 
2 A log C for quantifying how signihcantly a source emerges from the background, comparing 
the likelihood function C with and without that source. 


3.1. Detection and Localization 


This section describes the generation of a list of candidate sources, with locations and 
initial spectral hts, for processing by the standard LAT science analysis tools, espe cially 
gtlik e to compute the likelihood (§ 13.21) . This initial stage uses instead pointlike (IKerr 
20101) . Gompared with the gtlike-hased analysis described in 13.21 to ITTI it uses the same 


data, exposure, and IRFs, but the partitioning of the sky, the computation of the likelihood 
function, and its optimization, are independent. Since this version of the computation of 
the likelihood function is used for localization, it needs to represent a valid estimate of the 
probability of observing a point source with the assumed spectral function. 


The process started with an initial set of sources from the 2FGL analysis; not just those 
reported in that catalog, but also including all candidates failing the signihcance threshold 
(i.e., with TS < 25). It also used the latest extended source list with 25 entries (§ 13.4p . 
and the three-source representation of the Grab (§ 13.31) . The same spectral models were 
considered for each source as in § 13.31 but the favored model (power law or curved) was not 
necessarily the same. 


^See http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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Many details o f the processing were identical to the 2FGL procedure: using HEALPix@ 
flGorski et ahM2005[) with Ngide = 12, to tile the sky, resulting in 1728 tiles of ~25 deg^ area; 
optimizing spectral parameters for the sources within each tile, for the data in a cone of 
5° radius about the center of the tile; and including the contributions of all sources within 
10° of the center. The tiles are of course discrete, but the regions, which we refer to as 
Rols, for Regions of Interest, are overlapping and not independent. The data were binned 
in energy (14 energy bands from 100 MeV to 316 GeV) and position, where the angular bin 
size (the bins also defined using HEALPix) was set to be small compared with the PSF for 
each energy, and eyent type. Separating the photons according to eyent type is important, 
especially for localization, since the Front-converting events have a factor of two narrower 
PSF core than the Back-converting events. Thus the parameter optimization was performed 
by maximizing the logarithm of the likelihood, expressed as a sum over each energy band 
and each of the two event types {Front, Back). The hts for each Rol, maximizing the 
likelihood as a function of the free parameters, were performed independently. Gorrelations 
between sources in neighboring Rols were then accounted for by iterating all 1728 fits until 
the changes in the log likelihoods for all Rols were less than 10. 


After a set of iterations had converged, then the localization procedure was applied, 
and source positions updated for a new set of iterations. At this stage, new sources were 
occasionally added using the residual TS procedure described below. The detection and 
initial localization process resulted in 4029 candidate point sources with TS > 10. 


New features that are discussed below include an assessment of the reliability of each 
spectral fit and of the model as a whole in each Rol, a different approach to the normalization 
of the Galactic diffuse background component, and a method to unweight the likelihood to 
account for the effect of potential systematic errors in the Galactic diffuse on source spectra. 


3.1.1. Fit validation 

An important criterion is that the spectral and spatial models for all sources are not 
only optimized, but that the predictions of the models are consistent with the data. Maxi¬ 
mizing the likelihood as an estimator for spectral parameters and position is valid only if the 
likelihood, given a set of parameters, corresponds to the probability of observing the data. 
We have three measures. 

The hrst compares the number of counts in each energy band, combining Front and 


^http://healpix.sourceforge.net 
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Back, for each of the 1728 regions, dehning a x^-like measure as the sum of the squares of the 
deviations divided by the predicted number of counts. The number of counts is the expected 
variance for Poisson counting statistics. This measure is of course only a component of the 
likelihood, and depends only weakly on most of the point sources. That is, maximizing 
the likelihood does not necessarily minimize this quantity. But it is important to check 
the reliability of the diffuse model used, since this can distort the point source spectral hts. 
Figure [2] shows the distribution of that y^-like measure and its values as a function of location 
on the sky. The number of degrees of freedom is 14 (the number of energy bands) minus the 
effective number of variables. The fact that the distribution peaks at ~ 9 seems sensible. 
The ~35 regions with > 50 indicate problems with the model. Most are close to the 
Galactic plane, indicating difficulty with the component representing the Galactic diffuse 
emission. The few at high latitudes could be due to missing sources or, for very strong 
sources, inadequacy of the simple spectral models that we use. 



Fig. 2.— Distributions of the measure of consistency of the measured spectrum of each 
Rol with the model (capped at 50). Left: histogram highlighting the low-latitude subset. 
Right: distribution of the values over the sky. 

The second measure is a check that the spectral model for each source is consistent with 
the data. The likelihood associated with a source is the product of the likelihoods for that 
source for each energy band, including the contributions of nearby, overlapping sources, and 
the diffuse backgrounds. The correlations induced by those are only relevant for the lower 
energies, typically below 1 GeV. For this analysis, we keep these contributions hxed. We 
form the spectral £t quality as 2 log(£bands/'^fit) where the flux for each band is optimized 
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independently in £bands whereas the spectral model is applied in £fit. The spectrum in 
Figure [3] illustrates the concept. 

In Figure IU we show the distribution of the spectral £t quality for all preliminary 
spectra, with separate plots for the three different spectral functions (§133]): power law, log¬ 
normal, and power law with an exponential cutoff. The latter, applied almost exclusively to 
pulsars, is separated into sources in and out of the Galactic plane. It is seen that sources in 
the plane often have poorer hts. All are compared with an example distribution with 10 
degrees of freedom. There are 14 bands, and two to four parameters, but the higher-energy 
bands often do not contribute, so the number of degrees of freedom is not well defined and 
we use 10 for illustration only. 



Energy (GeV) 

Fig. 3.— The spectral energy distribution for a typical source, in this case PSR J1459—6053, 
as measured by the pointlike analysis. The lower plot shows the pulls, dehned as the square 
root of the difference 2 A log C between the htted flux and the spectral model in each energy 
band, signed with the residual. The points with error bars reflect the dependence of the 
likelihood on the flux for each energy band, combining Front and Back, while the curve is 
the result of the fit to all the energy bands. 

Finally, the localization process hts the logarithm of the likelihood as a function of 
position to a quadratic form, and checks the consistency with a y^-like measure (§ I3.1.3p . 
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Fig. 4.— Distributions of the spectral fit quality (capped at 50). Left: sources fit with a 
power-law spectrum; Center: sources £t with a log-normal; Right: sources fit with a power 
law with exponential cutoff. All are overlaid with the distribution with 10 degrees of 
freedom. 


3.1.2. Galactic diffuse normalization and unweighting 

The model that we used for the Galactic diffuse background is a global £t using the 
data, as described in § 12.31 For an individual Rol however, we found that we needed to 
adjust the normalizations for each band to £t the data. For the relatively broad energy 
bands, four per decade, used in the pointlike £t we allow the normalization for each band to 
vary, effectively ignoring the spectral prediction of the diffuse component analysis. So, for 
each of the 1728 Rols, and for each of the eight energy bands below 10 GeV, we measured a 
normalization factor, which applies to both Front and Back, by maximizing the likelihood 
with respect to it. A motivation for this procedure was that, for the lowest energy bands, it 
often improved the £t consistencies of the spectral models of the sources in the same Rol. 

While the precision of the determination of the average contribution from the Galac¬ 
tic diffuse for an energy band is subject to only the statistics of the number of photons, 
the value of the Galactic diffuse intensity at the location of each source, that is, the an¬ 
gular distribution of the intensity, is subject to an additional systematic error. Since this 
intensity is strongly correlated with the measurement of the flux from the source itself, and 
the correlation can be very significant for weak sources, we have adopted an ad hoc, but 
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conservative procedure to account for the additional uncertainty by increasing the width 
of the log likelihood distribution from each energy band according to how sensitive it is 
to the Galactic diffuse contribution. This is accomplished by dividing the log likelihood by 
max(l,A^"difr/1000) where N^m is the predicted number of Galactic diffuse photons in the Rol. 
This has the effect of limiting the precision to the statistics of 1000 photons in the Rol and 
energy band, i.e. it unweights contributions from energy ranges for which the contribution 
from the diffuse component is relatively less well dehned. 


3.1.3. Localization 

The position of each source was determined by maximizing the likelihood with respect 
to its position only. That is, all other parameters are kept hxed. The possibility that a 
shifted position would affect the spectral models or positions of nearby sources is accounted 
for by iteration. Ideally, the likelihood is the product of two Gaussians in two orthogonal 
angular variables. Thus the log likelihood is a quadratic form in any pair of angular variables, 
assuming small angles. We dehne LTS, for Localization Test Statistic, to be twice the log 
of the likelihood ratio of any position with respect to the maximum; the LTS evaluated for 
a grid of positions is called an LTS map. We £t the distribution of LTS to a quadratic 
form to determine the uncertainty ellipse, the major and minor axes and orientation. We 
also dehne a measure, the localization quality (LQ), of how well the actual LTS distribution 
matches this expectation by reporting the sum of the squares of the deviations of eight points 
evaluated from the £t at a circle of radius corresponding to twice the geometric mean of the 
two Gaussian sigmas. Figure [5] shows examples of localization regions for point sources. The 
distribution of the localization quality is shown in Figure El 

An important issue is how to treat apparently significant sources that do not have good 
localization fits, which we dehned as LQ > 8. An example is shown in Figure |5] (right). We 
hagged such sources (Flag 9 in Table [3]) and for them estimated the position and uncertainty 
by performing a moment analysis of the LTS function instead of htting a quadratic form. 
Some sources that did not have a well-dehned peak in the likelihood were discarded by hand, 
on the consideration that they were most likely related to residual dihuse emission. Another 
possibility is that two nearby sources produce a dumbbell-like shape; for some of these cases 
we added a new source by hand. A hnal selection demanding that the semi-major radius 
(la) be less than 0?25 resulted in 3976 candidate sources of which 142 were localized using 
the moment analysis. 

As in IFGL and 2FGL, we compared the localizations of the brightest sources with 
associations with their true positions in the sky. This indicated that the absolute precision 
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Fig. 5.— Examples of localization TS maps. The contours for 68%, 95%, and 99% con¬ 
tainment are shown. The scale (in decimal degrees) is not the same in both plots. Left: 
PSR J1459—6053, a good localization with LQ = 0.63. Right: 3FGL J2246.7—5205, a bad 
localization with LQ = 14. 



Fig. 6.— The distribution, in the preliminary source list, of the localization quality LQ 
(capped at 10). 
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is still the same, ~0?005 at the 95% conhdence level. After the associations procedure (l d5.2p . 
we compared the distribution of distances to the high-conhdence counterparts (in units of 
the estimated la errors) with a Rayleigh distribution, and noted that it was slightly broader, 
by a factor 1.05 (smaller than the 1.1 factor used in IFGL and 2FGL). Gonsequently, we 
multiplied all error estimates by 1.05 and added 0?005 in quadrature to both 95% ellipse 
The result i ng co mparison with the Rayleigh distribution is shown in Figure 3 of 


axes. 


Ackermann et al.l (120151. SLAG) and indicates good agreement. 


3.1.4- Detection of additional sources 

We used the pointlike definition of likelihood itself to detect sources that needed to be 
added to the model of the sky. Using HEALPix with Wide = 512, we defined 3.2 M pixels in 
the sky, separated by ~ 0?15, then evaluated the improvement in the likelihood from adding 
a new point source at the center of each, assuming a power-law spectrum with index 2.2. 
The TS value for each attempt, assigned to the pixel, defines a residual TS map of the sky. 
Next we performed a cluster analysis for all pixels with TS > 10, determining the number 
of pixels, the maximum TS, and the TF-weighted centroid. All such clusters with at least 
two pixels were added to a list of seeds. Then each seed was reanalyzed, now allowing the 
spectral index to vary, with a full optimization in the respective Rol, and then localized. 
The last step was to add all such refit seeds, if the hts to the spectrum and the position were 
successful, and TS > 10, as new sources, for a final optimization of the full sky. 


3.2. Significance and Thresholding 

The framework for this stage of the analysis is inherited from the 2FGL catalog. It 
splits the sky into Rols, varying typically half a dozen sources near the center of the Rol at 
the same time. There were 840 Rols for 3FGL, listed in the ROIs extension of the catalog 
(App. 1^. The global best fit is reached iteratively, injecting the spectra of sources in the 
outer parts of the Rol from the previous step. In that approach the diffuse emission model 
(§ 12.31) is taken from the global templates (including the spectrum, unlike what is done 
with pointlike in § 13.ip but it is modulated in each Rol by three parameters: normalization 
and small corrective slope of the Galactic component and normalization of the isotropic 
component. Appendix lAI shows how those parameters vary over the sky. 

Among more than 4000 seeds coming from the localization stage, we keep only sources 
at US' > 25, corresponding to a signihcance of just over 4cr evaluated from the distribution 
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with 4 degrees of freedom (position and spectral parameters. iMattox et al.lll996l) . The model 
for the current Rol is readjusted after removing each seed below threshold, so that the hnal 
model fits the full data. The low-energy flux of the seeds below threshold (a fraction of 
which are real sources) can be absorbed by neighboring sources closer than the PSF radius. 
There is no pair of seeds closer than 0?1, so the neighbors are unaffected at high energy. The 
fixed sources outside the core of the Rol are not tested and therefore not removed during 
the last fit of an Rol. Since the TS threshold at the previous step was set to 16, seeds with 
16 < TR < 25 still populate the outer parts of the Rol, preventing the background level to 
rise (bullet 5 below). 

We introduced a number of improvements with respect to 2FGL (by decreasing order 
of importance): 


1. After 2FGL was completed we understood that it was important to account for the 
different instrumental backgrounds in Front and Back events flj I2.3p . Implicitly as¬ 
suming that they were equal as in 2FGL resulted in lower TS (fewer sources) and 
tended to underestimate the low-energy flux. The impact is largest at high latitude. 
We used different isotropic spectral templates for Front and Back events, but a com¬ 
mon renormalization parameter. We also used different Front and Back models of the 
Earth limb. The same distinction was introduced for computing the fluxes per energy 
band (§ 13.5p and per month (§ 13.6p . 


2. Another effect discovered after 2FGL was a slight inconsistency (8% at 100 MeV) 
between the Front and Back effective areas. This affected mostly the Galactic plane, 
where the strong interstellar emission makes up 90% of the events. That effect created 
opposite low-energy residuals in Front and Back which did not compensate each other 
because of the differi ng PSF. It was corre cted empirically in the P7REP_SOURGE_V15 
version of the IRFs flBregeon et al.l 120131) . 


3. We put in place an automatic iteration procedure at the next-to-last step of the process 
checking that the all-sky result is stable (2FGL used a fixed number of five iterations), 
similar to what was done for localization in 2FGL. Quantitatively, we iterated an Rol 
and its neighbors until log£ did not change by more than 10. In practice this changes 
nothing at high latitude, but improves convergence in the Galactic plane. Fifteen 
iterations were required to reach full convergence. That iteration procedure was run 
twice, allowing sources to switch to a curved spectral shape (§ 13.3p after the hrst 
convergence. 


4. The software issue which prevented using unbinned likelihood in 2FGL was solved. We 
took advantage of that by using unbinned likelihood at high energy where keeping track 
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of the exact direction of each event helps. At low energy we used binned likelihood in 
order to cap the memory and CPU resources. The dividing energy was set to 3 GeV, 
resulting in data cubes (below 3 GeV) and event lists (above 3 GeV) of approximately 
equal size. Both data sets were split between Front and Back. This was implemented 
in the SummedLikelihood framework of pyLikelihood. In binned mode, the pixel size 
was set to 0?2 and 0?3 for Front and Back events, respectively (at 3 GeV the full width 
at half maximum of the PSF is 0?25 and 0?38, respectively). The energy binning was 
set to 10 bins per decade as in 2FGL. In the exposure maps for unbinned mode, the 
pixel size was set to 0?1 (even though the exposure varies very slowly, this is required 
to model precisely the edge of the held of view). 

5. We changed the criterion for including sources outside the Rol in the model. We 
replaced the hat 7° distance threshold by a threshold on contributed counts (predicted 
from the model at the previous step). We kept all sources contributing more than 2% of 
the counts per square degree in the Rol. This is a good compromise between reliability 
and memory/GPU requirements, and accounts for bright sources even far outside the 
Rol (at 100 MeV the 95% containment radius for Back events is 14°). Gompared to 
2FGL, that new procedure ahects mostly high latitudes (where the sources make up 
a larger fraction of the dihuse emission). Because it brings more low-energy events 
from outside in the model, it tends to reduce the htted level of the low-energy dihuse 
emission, resulting in slightly brighter and softer source spectra. 

6. The hts are now performed up to 300 GeV, and the overal signihcances (Signif_Avg) 
as well as the spectral parameters refer to the full 100 MeV to 300 GeV band. 

7. We introduced explicitly the model of the Sun and Moon contributions 12.3p . with¬ 
out any adjustment or free parameter in the likelihood analysis. The success of that 
procedure is illustrated in Figures 171 and ITOl 

8. For homogeneity (so that the result does not depend on which spectral model we start 
from) the TS > 25 threshold was always applied to the power-law model, even if the 
best-ht model was curved. There are 21 sources in 2FGL with TS — TScmve < 25 which 
would not have made it with this criterion (see li 13.31 for the dehnition of TScuwe)- 


3.3. Spectral Shapes 


The spectral representation of sources was mostly the same as in 2FGL. We introduced 
an additional p arameter modelin g a super- or subexponentially cutoff power law, as in the 
pulsar catalog flAbdo et al.l 120131) . However this was applied only to the brightest pulsars 
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(PSR J0835-4510 in Vela, J0633+1746, J1709-4429, J1836+5925, J0007+7303). The global 
fit with nearby sources was too unstable for the fainter ones, which were left with a simple 
exponentially cutoff power law. The subexponentially cutoff power law was also adopted 
for the brightest blazar 3C 454.The £t was very significantly better than with either a 
log-normal or a broken power law shape. Even though bright sources are not a scientihc 
objective of a catalog, avoiding low-energy spectral residuals (which translate into spatial 
residuals because of the broad PSF) is important for nearby sources. 

Therefore the spectral representations which can be found in 3FGL are: 


• a log-normal representation (LogParabola in the tables) for all signihcantly curved 
spectra except pulsars and 3C 454.3: 


diV 

dF 


K 


/ -a-P\og(E/Eo) 

\Eo) 


( 1 ) 


where log is the natural logarithm. The reference energy Eq is set to Pivot_Energy in 
the tables. The parameters K, a (spectral slope at Eq) and the curvature jS appear 
as Flux_Density, Spectral_Index and beta in the tables, respectively. No negative 
(3 (spectrum curved upwards) was found. The maximum allowed jS was set to 1 as in 
2FGL. 


• an exponentially cutoff power law for all signihcantly curved pulsars and a super- or 
subexponentially cutoff power law for the bright pulsars and 3G 454.3 (PLExpCutoff 
or PLSuperExpCutoff in the tables, depending on whether b was hxed to 1 or left free): 



where the reference energy Eo is set to Pivot_Energy in the tables and the parameters 
K, T (low-energy spectral slope), Ec (cutoff energy) and b (exponential index) appear 
as Flux_Density, Spectral_Index, Cutoff and Exp_Index in the tables, respectively. 
Note that this is not the way that spectral shape appears in the Science Tools (no 
{Eq/Ec)^ term in the exponential), so the error on K in the tables was obtained from 
the covariance matrix. The minimum T was set to 0.5 (in 2FGL it was set to 0). 

• a simple power-law form for all sources not significantly curved. 


®That is only a mathematical model, it should not be interpreted in a physical sense since it is an average 
over many different states of that very variable object. 
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As in 2FGL, a source is considered significantly curved if TS'curve > 16 where TScmve = 
2(log£(curved spectrum)—log £ (power-law)). The curved spectrum is PLExpCutoff (or 
PLSuperExpCutoff) for pulsars and 3C 454.3, LogParabola for all other sources. The cur¬ 
vature significance is reported as Signif_Curve (see 13.5p . 

Another difference with 2FGL is that the complex spectrum of the Grab was represented 
as three components: 


a PLExpCutoff shape for the pulsar, with free K, F and Ec¬ 


a soft power-law shape for the synchrotron en iission of the nebu la, with free K and 
F since the synchrotron emission is variable flAbdo et al.l l2011a) . The synchrotron 
component is called 3FGL J0534.5-|-2201s. 

a hard power-law shape for the inverse Gompto n emission of the nebula, with param¬ 
eters fixed to those found in lAbdo et al.l fl2010el) . That component does not vary, and 
leaving it free made the fit unstable. It is called 3FGL J0534.5-|-2201i. 


In 2FGL, two sources (MSH 15—52 and Vela X) spatially coincident with pulsars had 
trouble converging and their spectra were fixed to the result of the dedicated analysis 
flAbdo et al.N2010alJd) . In 3FGL the spectra of five sources were fixed for the same rea¬ 
son: th e same two, the Inverse Gompton component of the Grab Nebula, the Gygnus X 
cocoon flAckermann et ahll2011al) and th e y-Gy gni supernova remnant. The spatial template 
of y-Gygni was taken from Lande et al. f 2ni3l as in IFHL. We did not switch to the more 
complex spatial template used in I Ackerman n et ahl fcOllal) but the spectral template was 
obtained from a reanalysis of the Gygnus region including the Gygnus X cocoon (L. Tibaldo, 
private communication). 


Overall in 3FGL six sources (the five brightest pulsars and 3G 454.3) were fit as 
PLSuperExpCutoff (with b of Eq. [2] < 1), 110 pulsars were fit as PLExpCutoff, 395 sources 
were fit as LogParabola and the rest (including the five fixed sources) were represented as 
power laws. 


3.4. Extended Sources 

As for the 2FGL and IFHL catalogs, we explicitly model as spatially extended those 
LAT sources that have been shown in dedicated analyses to be resolved by the LAT. Twelve 
extended sources were entered in the 2FGL catalog. That number grew to 22 in the IFHL 
catalog. The spatial templates were based on dedicated analysis of each source region and 
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have been normalized to contain the entire flnx from the sonrce (> 99% of the flnx for 
nnlimited spatial distribntions snch as 2-D Ganssians). The spectral form chosen for each 
sonrce is the best adapted among those nsed in the catalog analysis (see I3.3p . Three more 
extended sonrces have been reported since then and were inclnded in the same way in the 
3FGL analvsiJ^. 


The catalog process does not involve looking for new extended sonrces or testing pos- 
sible extensi o n of sonrces detected as point-like. This was last done comprehensively by 


Lande et ahl 020121 ) based on IFGL. The extended sonrces published since then were the 


result of focussed studies so there most likely remain unreported faint extended sources in 
the Fermi-LAT data set. The process does not attempt to reht the spatial shape of known 
extended sources either. 


The extended sources include twelve supernova remnants (SNRs), nine pulsar wind 
nebulae (PWNe) or candidates, the Gygnus X cocoon, the Large and Small Magellanic 
Glouds (LMG and SMG), and the lobes of the radio galaxy Gentaurus A. Below we provide 
notes on new sources and changes since 2FGL: 


HB 21 is an SNR recently reported as a LAT source (IReichardt et al.ll2012l) . We added 
it to the list, us i ng the simple disk template and LogParabola spectral shape derived 
by Pivato et ah ( 2013 ). 


HESS J1303—631 and HESS J1841—055 are two H.E.S.S. sources (m ost likely PWNe) 
recently reported as faint hard LAT sources bv lAcero et aLl (1201311 . We added them 
to the list, using the original H.E.S.S. template rather than the best spatial ht to the 
LAT data, in keeping with the spectral analysis in that paper. 


• We changed the spectral representation of the LMG and the Gygnus Loop from PLExpCutof f 
to LogParabola, which hts the data better. The curvature of the fainter SMG spectrum 
is not signihcant; therefore it was ht as a power law. 


In general, we did not allow any point source inside the extended templates, even when 
the TS maps indicated that adding new seeds would improve the ht. Most likely (pending 
a dedicated reanalysis) those additional seeds were simply residuals due to the fact that the 
very simple geometrical representations that we adopted are not precise enough, rather than 
independent point sources. We preferred not splitting the source hux into pieces. The only 
exceptions are 3FGL J1823.2—1339 within HESS J1825—137, 3FGL J2053.9-1-2922 inside 


^°The templates and spectral models are available through the Fermi Science Support Center. 
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the Cygnus Loop, 3FGL J0524.5—6937 inside the LMC, and sonrces inside the Cygnus X 
cocoon. The hrst one is as signihcant as the extended sonrce and was a 2FGL source already. 
The next two are well localized over large extended sources and show a very hard spectrum, 
so they do not impact the spectral characteristics of the extended sources. The Gygnus X 
cocoon was fixed (§ 13.3p and allowing point sources on top of it was necessary to reach a 
reasonable representation of the region. 

Table [T] lists the source name, spatial template description, spectral form and the ref¬ 
erence for the dedicated analysis. These sources are tabulated with the point sources, with 
the only distinction being that no position uncertainties are reported and their names end 
in e (see § 14.ip . Unidentified point sources inside extended ones are marked by “xxx held” 
in the ASS0C2 column of the catalog. 


3.5. Flux Determination 

The source photon huxes are reported in the same hve energy bands (100 to 300 MeV; 
300 MeV to 1 GeV; 1 to 3 GeV; 3 to 10 GeV; 10 to 100 GeV) as in 2FGL. The huxes 
were obtained by freezing the spectral index to that obtained in the ht over the full range 
and adjusting the normalization in each spectral band. For the curved spectra (§ 13.3p the 
spectral index in a band was set to the local spectral slope at the logarithmic mid-point 
of the band y/EnEn+i, restricted to be in the interval [0,5]. The photon hux between 1 
and 100 GeV as well as the energy hux between 100 MeV and 100 GeV (F 35 and S25 in 
Table [SI the subscript ij indicates the energy range as 10*-10-^ MeV), are derived from 
the full-band analysis assuming the best spectral shape, and their uncertainties from the 
covariance matrix. Even though the full analysis is carried out up to 300 GeV in 3FGL, we 
have not changed the energy range over which we quote huxes so that they can be easily 
compared with past huxes. The photon hux above 100 GeV is negligible anyway and the 
energy hux above 100 GeV is not precisely measured (even for hard sources). 

Improvements with respect to the 2FGL analysis are: 

• We used binned likelihood in the hrst three bands (up to 3 GeV) and unbinned likeli¬ 
hood in the last two bands, distinguishing Eront and Back events. The pixel sizes in 
each band in binned mode were 0?3 and 0?5, 0?2 and 0?3, 0?1 and 0?15 where in each 
band, the hrst value is for Eront, the second one for Back. This reduces error bars by 
10-15% compared to mixing Eront and Back events as in 2FGL. 

• Following what was done in the IFHL catalog, the errors on the huxes of moderately 
faint sources (TF > 1 in the band) were computed as la errors with MINOS in the 
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Table 1. Extended Sources Modeled in the 3FGL Analysis 


3FGL Name 

Extended Source 

Spatial Form 

Extent (deg) 

Spectral Form 

J0059.0-7242e 

SMC 

2D Gaussian 

0.9 

PowerLaw 

J0526.6-6825e 

LMC 

2D Gaussian^ 

1.2, 0.2 

LogParabola 

J0540.3+2756e 

S 147 

Map 


PowerLaw 

J0617.2+2234e 

IC 443 

2D Gaussian 

0.26 

LogParabola 

J0822.6-4250e 

Puppis A 

Disk 

0.37 

PowerLaw 

J0833.1-4511e 

Vela X 

Disk 

0.88 

PowerLaw 

J0852.7-4631e 

Vela Junior 

Disk 

1.12 

PowerLaw 

J1303.0-6312e 

HESS J1303-631 

2D Gaussian 

0.16 

PowerLaw 

J1324.0-4330e 

Centaurus A (lobes) 

Map 


PowerLaw 

J1514.0-5915e 

MSH 15-52 

Disk 

0.25 

PowerLaw 

J1615.3-5146e 

HESS J1614-518 

Disk 

0.42 

PowerLaw 

J1616.2-5054e 

HESS J1616-508 

Disk 

0.32 

PowerLaw 

J1633.0-4746e 

HESS J1632-478 

Disk 

0.35 

PowerLaw 

J1713.5-3945e 

RX J1713.7-3946 

Map 


PowerLaw 

J1801.3-2326e 

W28 

Disk 

0.39 

LogParabola 

J1805.6-2136e 

W30 

Disk 

0.37 

LogParabola 

J1824.5-1351e 

HESS J1825-137 

2D Gaussian 

0.56 

LogParabola 

J1836.5-0655e 

HESS J1837-069 

Disk 

0.33 

PowerLaw 

J1840.9-0532e 

HESS J1841-055 

2D Gaussian^ 

(0.41, 0.25) 

PowerLaw 

J1855.9+0121e 

W44 

Ring*^ 

(0.22, 0.14), (0.30, 0.19) 

LogParabola 

J1923.2+1408e 

W51C 

Disk*^ 

(0.40, 0.25) 

LogParabola 

J2021.0+4031e 

7 -Cygni 

Disk 

0.63 

PowerLaw 

J2028.6+4110e 

Cygnus X cocoon 

2D Gaussian 

2.0 

PowerLaw 

J2045.2+5026e 

HB 21 

Disk 

1.19 

LogParabola 

J2051.0+3040e 

Cygnus Loop 

Ring 

0.7, 1.6 

LogParabola 


Reference 


I<atsut^^^aL 

^and^^^aL ( 2 ^\^ 

^^nak^^^aL ) 
Aharoniar^^aL ^005) 
(2^^c) 

^and^^^aL ( 2 ^\^ 
^and^^^aL ( 2 ^\^ 
^and^^^aL ( 2 ^\^ 

^ell^^^al^ (2^^) 

^rondii^^aL ) 

^and^^^aL ( 2 ^\^ 
^laroniai^^aL ^ 008 ) 

^009b) 
^and^^^aL ( 2 ^\^ 
Ackermam^^aL 

I^atagir^^al^ 


^Combination of two 2D Gaussian spatial templates. 

^The shape is elliptical; each pair of parameters (a, b) represents the semi-major (a) and semi-minor (6) axes. 


Note. — List of all sources that have been modeled as extended sources. The Extent column indicates the radius for Disk sources, 
the dispersion for Gaussian sources, and the inner and outer radii for Ring sources. 
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Fig. 7.— Spectral energy distribution of the same source (the BL Lac TXS 0836+182) 
as 2FGL J0839.4+1802 (left) and 3FGL J0839.6+1803 (right). This source was flagged as 
‘possibly due to the Sun’ (Flag 11 in Table [3]) in 2FGL. Entering the Sun and Moon into 
the background model has reduced to zero the low energy signal that drove the 2FGL fit, 
resulting in a hard source in 3FGL. The dashed line is the best fit over the full energy range 
fl d3.2p and the gray-shaded area is the statistical uncertainty around the best fit (for a given 
spectral form). The vertical scale is not the same in the left and right plots. Note that the 
300 MeV to 1 GeV point has (asymmetric) error bars in 3FGL as explained in § 13.51 even 
though its significance is less than 2 ct. Upper limits (indicated by a downward triangle in 
3FGL and a downward arrow in 2FGL) are at 95% confidence level. 
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Energy (GeV) 

Fig. 8.— Spectral energy distribution of 3FGL J0222.6+4301 (3C 66A) fitted by a 
LogParabola spectrum with = 0.039 + 0.007 but Signif_Curve (defined just before Eq. [3]) 
= 2.81. The curvature is statistically significant but a power law cannot be excluded given 
the range of the systematic errors on effective area. 
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Energy (GeV) 

Fig. 9.— Spectral energy distribution of 3FGL J1226.9—1329 (PMN J1226—1328) flagged 
with bad spectral fit quality (Flag 10 in Table [3]). The hrst two points deviate from the 
power-law £t. This source is within 1?3 of the much brighter pulsar PSR J1231—1411 so it 
is confused with it (within rgg) below 600 MeV. 
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Minuit0 package. This was done whenever the relative error on flnx in the qnadratic 
approximation (from the covariance matrix) was larger than 10%. Both errors (lower 
and npper) are reported in the FITS table (App.[B|). The lower error is reported with a 
minus sign (when the error comes from the quadratic approximation, the lower error is 
simply minus the upper error). The upper limits for very faint so urces {TS < 1) 
were computed as in 2FGL, using the Bayesian method flHelenel 1198311 at 95% of the 
posterior probability. The upper error is then reported as 0.5{F^^ — F^^) where F^^ 
is the best-£t flux, and the lower error is set to NULL. 


• The same iteration procedure described in § 13.21 was put in place for the fluxes per 
energy band using a more stringent criterion (Alog£ < 3). Convergence was fast at 
high energy (little cross-talk between sources). It was a little slower at low energy (6 
iterations in the hrst band) but much faster than the full-band £t because no spectral 
adjustment was involved. 


• We report as nuFnuxxx_yyy the Spectral Energy Distribution (SED) in the band defined 
by XXX to yyy MeV, which can be directly overlaid on an SED plot. The SED was 
obtained by dividing the energy flux in the band by the band width in natural logarithm 
log(yyy/xxx). Since the £t is performed on the flux only (no spectral freedom in each 
band), the relative error on the SED is the same as that on the corresponding flux. 


As in 2FGL we report in 3FGL a curvature signihcance Signif_Curve = 

(in a units) after approximately accounting for systematic uncertainties on effective area via 

p _ E.(f. - 

where i runs over all bands, Ff^ is the flux predicted by the power-law model and F®* is the 
flux predicted by the best-fit (curved) model in that band from the spectral fit to the full 
band, reflects the systematic uncertainty on effective area (§ 13.7p . The values were set 
to 0.1, 0.05, 0.05, 0.05, 0.1 in our five bands (the fourth one went down from 0.08 in 2FGL, 
thanks to improved calibration). Eq. [3] is not exactly the same formula used for 2FGL. In 
2FGL Ff^ would have been replaced by Ff^. The disadvantage of the previous estimate was 
that it capped Signif_Curve to rather low values (below 15) resulting in a small dynamic 
range because the largest relative systematic errors are in the two extreme bands and in 
those bands the power-law fit can run way above the points (because the spectra are curved 
downwards). Using the curved fit (closer to the points) to estimate the systematic errors is 
a more reasonable procedure, and recovers a larger dynamic range (up to 85 in 3FGL). 


^^http://Icgapp.cern.ch/project/cls/work-packages/mathlibs/minuit/home.html 
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As in 2FGL we consider that only sources with Signif_Curve > 4 are signihcantly 
curved (at the 4 ct level). When i?syst is small (bright source) it can happen that TS'curve > 16 
(triggering a curved model following § I3.3p but Signif_Curve < 4. The 43 such sources 
with LogParabola spectra (and 2 pulsars with PLExpCutoff spectra) but Signif_Curve < 
4 could be power laws within systematic errors. Nevertheless we do not go back to power- 
law spectra for those sources because they are better fit with curved models and power-law 
models would result in negative low-energy residuals which might affect nearby sources. One 
of them is illustrated in Figure |H1 All are bright sources with modest curvature. 


Spectral £t quality (for Flag 10 in Table [3]) is computed as in Eq. 3 of iNolan et ah 


fj2012l. 2FGL) rather than as in § 13.1.11 Among the 42 sources flagged because of a too large 
spectral fit quality, most show deviations at low energy and are in confused regions or close 
to a brighter neighbor, as in Figure O 


Spectral plots for all 3FGL sources overlaying the best model on the individual SED 
points are available from the FSSG. 


3.6. Variability 

The light curves were computed over the same (1-month) intervals as in IFGL and 
2FGL (there are now 48 points). The first 23 intervals correspond exactly to 2FGL. The 
fluxes in each bin were obtained by freezing the spectral parameters to those obtained in the 
fit over the full range and adjusting the normalization. We used unbinned likelihood over 
the full energy range for the light curves. Over short intervals it does not incur a large GPU 
or memory penalty and it preserves the full information. We used a different isotropic and 
Earth limb model for Front and Back events, as in the main fit (§ 13.2|) . We also used a 
different Sun/Moon model for each month (the Sun is obviously at a different place in the 
sky each month). That improvement, together with our removing the solar flares, effectively 
mitigated the peaks that we noted in the 2FGL light curves due to the Sun passage near the 
source (Flag 11 in Table [3]). We have not noted any obvious Sun-related peak in the 3FGL 
light curves (Figure [10]). 

As in the band fluxes calculation (§ 13.51) the errors on the monthly fluxes of moder¬ 
ately faint sources {TS > 1) were computed as lower and upper Icr errors with MINOS in 
Minuit. Both errors (lower and upper) are reported in the FITS table (Table [T3D so the 
Unc_Flux_History column is a 2 x 48 array. This allowed providing more information in the 
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MJD (days) 

Fig. 10.— Light curve of 3FGL J1315.7—0732 (NVSS J131552—073301) in the ecliptic plane. 
That source is signihcantly variable. The flares do not correspond to the times when the 
Sun passed through the region (vertical orange bands). The only effect of the Sun passage 
is somewhat larger error bars. The gray-shaded horizontal area materializes the systematic 
uncertainty of 2%. Upper limits (indicated by a downward triangle) are at 95% conhdence 
level. 
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Fig. 11.— Light curve of 3FGL J1616.2—5054e (HESS J1616—508). That is an extended 
source that should not be variable. Indeed the monthly fluxes are compatible with a constant, 
but not with the flux extracted over the full four years (dashed line with gray-shaded uncer¬ 
tainty). That inconsistency is due to a remaining difference between binned and unbinned 
likelihood fits affecting only extended sources. 



Flux (10 ph cm 



MJD (days) 

Fig. 12.— Light curve of 3FGL J2021.5+4026 (PSR J2021+4026 in the 7 Cygni SNR). The 
variability of that pulsar is easily detected by the automatic procedure. The vertical scale 
does not start at 0 . 
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light curve ploto by keeping points with error bars whenever TS > 1 (the lower error does 
not reach 0). When TS < 1 the 95% upper limit is converted into an upper error in the 
same way as in 2FGL and the band fluxes calculation. 

We noted an inconsistency between the light curve and the flux from the main fit (over 
the full interval) in several extended sources, whereby the average of the light curve appears 
distinctly above the flux from the main fit. It is particularly obvious in Cen A lobes, HESS 
J1616—508 fFigure ITT]) . S 147, W28, and W30. We traced the problem to the fact that 
we used unbinned likelihood over the whole energy range for the light curves, but binned 
likelihood for the main fit below 3 GeV. We have not found any evidence that this affects 
the point sources. Since we do not expect variability in extended sources, we left this 
inconsistency in the catalog as a known feature. 

The variability indicator Variability_Index is the same as in 2FGL, with the same 
relative systematic error of 2%. Variability is considered probable when Variability_Index 
exceeds the threshold of 72.44 corresponding to 99% confidence in a distribution with 47 
degrees of freedom. 


The Grab nebula and pulsar are a particularly difficult case. The nebula is very variable 
flTavani et al.l l201ll: lAbdo et al.l l201ld) while the pulsar has no detected variability. So we 
would have liked the synchrotron component to absorb the full variability in 3FGL. It does 
not turn out this way, however, because the spectrum of the nebula becomes much harder 
during flares. This is not accounted for in the variability analysis (the spectral slopes are 
fixed to that in the full interval). As a result, the pulsar component also increases during 
the nebular flares and the pulsar becomes formally variable. We stress here that it is only a 
feature of our automatic analysis and is in no way a real detection of variability in the Grab 
pulsar. Besides the Grab, we detect the (real) variability of PSR J2021+4026 (Figure YTR 


Allafort et al.ll2013[l . The only other formally variable pulsar is PSR J1732—3131 just above 


threshold. Since this is one in 137 pulsars, it is compatible with a chance occurrence at the 
99% confidence level. 
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Angular Distance (deg) 


Fig. 13.— Distribution of nearest-neighb or distances for 3FG L sources at |6| > 10°. The 
dashed curve was derived as described in lAbdo et all (12010(1 . IFGL). It is the distribution 
expected if sources could be detected at arbitrarily small angular separations. The dashed 
curve is normalized to match the total observed number of sources for separations > 0?8 
(2035). This corresponds to an expected true number of sources (extrapolated down to 0 
separation) of 2336 at |6| > 10°. 


3.7. Limitations and Systematic Uncertainties 

3.7.1. Source confusion 


As for the IFGL and 2FGL catalogs we investigated source confusion by comparing 
the actual distribution of angular separations between 3FGL sources with what would be 
expected for a population of sources that could be detected i ndependently regarcl less how 
small their angular separations. The formalism is dehned in lAbdo et al.l fl2010dl. IFGL). 


^^These plots are available from the FSSC. 
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We considered the region of the sky above |6| = 10°, within which the average angular 
separation of 3FGL sources is 2?2. The distribution of nearest-neighbor distances is shown 
in Figure Uni along with the distribution expected if the source detection efficiency did not 
decrease for closely-spaced sources. The observed density of nearest-neighbor starts to fall 
below the expected curve at about 0?8 angular separation. The implied number of missing 
closely-spaced sources is ~140, or about 6% of the estimated true source count in the region. 
For the 2FGL catalog the fraction was only 3.3%. This indicates that even though the PSF 
improved after the Pass 7 reprocessing, the larger number of detected sources (2193 vs. 1319) 
is now pushing the LAT catalog into the confusion limit even outside the Galactic plane. 
Because the confusion process goes as the square of the source density, the expected number 
of sources above the detection threshold within 0?5 of another one (most of which are not 
resolved) has increased by a factor of 3 between 2FGL and 3FGL. 

The consequence of source confusion is not only losing a fraction of sources. It can 
also lead to “composite” y-ray sources merging the characteristics of two very nearby astro¬ 
nomical objects. An example is the unassociated 3FGL J0536.4—3347, located between two 
bright blazars. Its spectrum is relatively soft, similar to that expected from the FSRQ BZQ 
J0536—3401, 14' away. Its location, however, is closer (4') to the BL Lac BZB J0536—3343 
because that one dominates at high energy where the Fermi PSF is best. That issue is 
discussed in more detail in the 3LAG paper. 


3.7.2. Instrument response functions 


The systematic uncertainties on effective area have improved since 2FGL, going from 
P7SOURGE_V6 to P7REP_SOURGE_V15. They are now estimated to be 5% between 
316 MeV and 10 GeV, increasing to 10% at 100 M eV and 15% at 1 T e V (see the caveats page 
at the FSSG), following the methods described bv I Acker maun et ahl fj2012al) . As in previous 
LAT catalogs, we have not included those uncertainties in any of the error columns, because 
they apply uniformly to all sources. They must be kept in mind when working with abso¬ 
lute numbers, but comparisons between sources can be carried out at better precision. The 
systematic uncertainties on effective area have been included in the curvature signihcance 
(§ 13.Sp and a systematic uncertainty of 2% on the stability of monthly flux measurements 
(measured directly on the bright pulsars) has been included in the variability index (§ 13.6p . 
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3.1.3. Diffuse emission model 


The model of diffuse emission is the main source of uncertainties for faint sources. 
Contrary to the effective area, it does not affect all sources equally: its effects are smaller 
outside the Galactic plane where the diffuse emission is fainter and varying on larger angular 
scales. It is also less of a concern in the high-energy bands (> 3 GeV) where the core of 
the PSF is narrow enough that the sources dominate the background under the PSF. But 
it is a serious concern inside the Galactic plane in the low-energy bands (< 1 GeV) and 
particularly inside the Galactic ridge (|/| < 60°) where the diffuse emission is strongest and 
very structured, following the molecular cloud distribution. It is not easy to assess precisely 
how large the uncertainties are, because they relate to uncertainties in the distributions of 
interstellar gas, the interstellar radiation held, and cosmic rays, which depend in detail on 
position on the sky. 


For an assessment we have tried re-extracting t he source spectr a using one of the eight 
alternative interstellar emission models described in Ide Palma et al.l (120131) . namely the one 
obtained with optically thin Hi, an SNR cosmic-ray source distribution and a 4 kpc halo, 
adapted to the P7REP IRFs. For computational reasons we have not used all eight alterna¬ 
tive models, but that one should be representative. In each Rol we left free the normalization 
of each component of the model contributing (with its normalization set to 1) more than 
3% of the total counts in the Rol. The isotropic normalization was also left free, but the 
inverse Gompton, Loop I and Fermi bubbles components were hxed (too large scale to be 
htted inside a single Rol). That approach (independent components) differs enough from the 
standard diffuse model that it can provide a stronger test than comparing with the previous 
generation diffuse model, as we did for 2FGL. Nevertheless both models still rely on nearly 
the same set of H l and GO maps of the gas in the interstellar medium, so they are not as 
independent as we would like. 


The results show that the systematic uncertainty more or less follows the statistical one, 
i.e., it is larger for fainter sources in relative terms. We list the induced biases and scatters 
of flux and spectral index in Table [2J We have not increased the flux and index errors in 
the catalog itself accordingly because this alternative model does not ht the data as well as 
the newer one. The ht quality is nearly everywhere worse, except near the Garina region 
where we know the standard model does not ht the data very well (App. |^. From that point 
of view we may expect these estimates of the systematic uncertainties to be upper limits. 
So we regard the values as qualitative estimates. In the Galactic plane (and even worse in 
the Galactic ridge) the systematic uncertainties coming from the dihuse model are larger 
than the statistical ones. In the Galactic ridge, even the bias is larger than the statistical 
uncertainty. The ehect is larger than what we estimated for 2FGL (even though the dihuse 
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model has improved), partly because the exposure is twice as deep and partly because the 
new alternative model is further from the standard one. Outside the Galactic plane the 
systematic uncertainty due to the diffuse model remains less than the statistical one, and 
the bias is negligible. 

The same comparison also allows flagging outliers as suspect f G3.9p . 119 sources received 
Flag 1 (Table [3]) because they ended up with TS <25 with the alternative model, and 118 
received Flag 3, indicating that their photon or energy fluxes changed by more than 3cr. That 
uncertainty also appears in Flag 4 whereby we flag all sources with source-to-background 
ratio less than 10% in all bands in which they are statistically signihcant. 


3.7.4- Analysis method 

The check presented in this section is new to 3FGL. As explained in 13. li the pointlike- 
based method used to detect and localize sources also provides an estimate of the source 
spectra and signihcance. Therefore we use it to estimate systematic errors due to the analysis 
itself. Many aspects differ between the two methods: the code, the Rols, the Earth limb 
representation. The alternative method does not remove faint sources (with TS < 25) from 
the model. The diffuse model is the same spatially but it was rescaled spectrally in each 
energy bin. The pointlike-hased method also rescales log C in order to play down the energy 
bins in which the source-to-background ratio is low. 

The procedure to compare the results is the same as in § 13.7.31 We list the induced 
biases and scatters of flux and spectral index in Table [2J In general, the effect of changing 
the analysis procedure is less than changing the diffuse model. Outside the Galactic ridge 
(and even outside the Galactic plane), we observe a negative bias on flux and index (i.e. 
fainter harder sources with the pointlike pipeline) close to half the statistical error. That 
effect is probably the result of removing the sources below threshold in the standard method. 
This favors absorbing the flux of faint neighbors at low energy (where the PSF is broad), 
resulting in somewhat brighter and softer sources. 

A total of 118 sources received Flag 1 (TS' < 25 with pointlike), and 101 received Flag 
3 (flux changed by more than 3 ct). Only 25 (Flag 1) and 19 (Flag 3) sources are flagged from 
both the diffuse model and the analysis method comparisons. In other words, the 3FGL 
catalog is more or less half way between the result from pointlike and the result with the 
alternative diffuse model. Gomparing the lists from pointlike and the alternative diffuse 
model would result in 202 sources with Flag 1 and 209 with Flag 3. 
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3.8. Sources Toward Local Interstellar Clouds and the Galactic Ridge 

As we did for the 2FGL catalog, we carefully evaluated which sources are potentially 
artifacts due to systematic uncertainties in modeling the Galactic diffuse emission. The pro¬ 
cedure, described in more detail in the 2FGL paper, flags unassociated sources with moderate 
TS and spectral index F > 2, corresponding to features in individual gas components. For 
3FGL we did not consider sources that have very curved spectra to be artifacts. Very soft 
sources with power-law spectra are instead more likely to be problematic. Sources considered 
to be potential artifacts are assigned an analysis flag in the catalog 13.9p . We also append 
c to the source names. 

Relative to the 2FGL catalog, far fewer c sources are flagged here (78 here vs. 162 
for 2FGL) despite the much greater number of sources overall in the 3FGL catalog. Away 
from the Galactic plane, the reduction of c sources is primarily due to improvement of 
the representation of the dark gas component of the Galactic diffuse emission model in the 
vicinity of massive star-forming regions 12 .3p . At low latitudes, the reduction primarily is 
due to relaxing the criterion on unassociated sources with very curved spectra. 

Figure [m shows the locations of the c sources for 3FGL. The majority are close to the 
Galactic plane, where the diffuse y-ray emission is brightest and very structured. Glusters are 
apparent in regions where spiral arms of the Milky Way are viewed essentially tangentially, 
in particular the Gygnus {I ~ 80°) and Garina {I ~ 285°) regions where the systematic 
uncertainties of the Galactic diffuse emission model are especially large. None of the c 
sources is identihed (§ 15.ip and 63 (~80%) have no hrm association with a counterpart 
at other wavelengths, a much larger fraction than the overall average (~30%) for 3FGL 
(Table E]). 


3.9. Analysis Flags 

As in 2FGL we identihed a number of conditions that should be considered cautionary 
regarding the reality of a source or the magnitude of the systematic uncertainties of its 
measured properties. They are described in Table Ei 

Each hag has the same dehnition as for the 2FGL catalog, except for Flag 7, which was 
unused in that catalog. 

Flags 1 to 12 have similar intent as in 2FGL, but diher in detail: 


Flags 1 and 3 are now applied not only when a source is sensitive to changing the 
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Selection 

Quantity 

Diffuse mode 

Bias 

i (§pmD 

Scatter 

Analysis methc 
Bias 

^d (§[X7 :td 

Scatter 

Galactic 

Ridge 

Eflux (174) 
Index (88) 

+1.88CT (+21%) 

+1.44ct (+0.14) 

3.40a (42%) 
1.81a (0.37) 

-0.47a (-7%) 
-0.08a (-0.01) 

1.93a (27%) 
2.40a (0.21) 

Galactic 

Plane 

Eflux (662) 
Index (470) 

+0.51a (+7%) 
+0.34a (+0.04) 

2.19a (32%) 
1.54a (0.21) 

-0.66a (-12%) 
-0.44a (-0.06) 

1.26a (23%) 
1.15a (0.15) 

High 

Latitude 

Eflux (2193) 
Index (1960) 

+0.07a (+1%) 
+0.23a (+0.03) 

0.98a (15%) 
0.73a (0.10) 

-0.42a (-7%) 
-0.34a (-0.05) 

0.74a (13%) 
0.73a (0.10) 


Table 2: The table gives the bias and the scatter induced by changing one of two important 
elements in the analysis chain, first in units of the statistical error (i.e., on {Af^ — Ai)/ai), 
then in absolute terms (i.e., on Af^ — Ai), where Ai is either the log of the energy flux 
between 100 MeV and 100 GeV or the spectral index in the standard analysis, Af^ is the 
same quantity in the alternative analysis and ai the statistical uncertainty on Ai. The 
spectral index comparison is restricted to pure power-law sources. The Galactic Ridge is 
defined as |6| < 2° and |/| < 60°. High Latitude is defined as \b\ > 10°. The Galactic Plane 
is everything else (i.e., it does not include the Galactic Ridge). The number of sources in 
each selection is given in parentheses after the quantity. 



Galactic Longitude 


Fig. 14.— Locations of the c sources in the 3FGL catalog overlaid on a grayscale represen¬ 
tation of the model for the Galactic diffuse y-ray emission used for the 3FGL analysis (see 
12.3p . The plotted symbols are centered on the locations of the sources. The model diffuse 
intensity is shown for 1 GeV and the spacing of the levels is logarithmic from 1% to 100% 
of the peak intensity. 
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Table 3. Definitions of the Analysis Flags 


Flag®^ Meaning 

1 Sonrce with TS > 35 which went to TS' < 25 when changing the diffuse model 

I3.7.3P or the analysis method I3.7.4p . Sources with TS <35 are not flagged 
with this bit because normal statistical fluctuations can push them to TS < 25. 

2 Not used. 

3 Flux (> 1 GeV) or energy flux (> 100 MeV) changed by more than 3 ct when 
changing the diffuse model or the analysis method. Requires also that the flux 
change by more than 35% (to not flag strong sources). 

4 Source-to-background ratio less than 10% in highest band in which TS > 25. 
Background is integrated over vrrgg or 1 square degree, whichever is smaller. 

5 Closer than from a brighter neighbor. 6'ref is defined in the highest band in 
which source TS > 25, or the band with highest TS* if all are < 25. 6'ref is set 
to 2? 17 (FWHM) below 300 MeV, 1?38 between 300 MeV and 1 GeV, 0?87 
between 1 GeV and 3 GeV, 0?67 between 3 and 10 GeV and 0?45 above 

10 GeV (2res). 

6 On top of an interstellar gas clump or small-scale defect in the model of 
diffuse emission; equivalent to the c designator in the source name flj 13.81) . 

7 Unstable position determination; result from gtfindsrc outside the 95% ellipse 
from pointlike. 

8 Not used. 

9 Localization Quality > 8 in pointlike (Ij 13.ip or long axis of 95% ellipse > 0?25. 

10 Spectral Fit Quality > 16.3 (Eq. 3 of Nolan et ah 2012. 2FGL). 

11 Possibly due to the Sun flj 13.61) . 

12 Highly curved spectrum; LogParabola f3 fixed to 1 or PLExpCutoff 
Spectral_Index fixed to 0.5 (see 13.3p . 


^In the FITS version the values are encoded as individual bits in a single column, with 
Flag n having value 2^”“^b For information about the FITS version of the table see Table ITBl 
in AppJHl 
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diffuse model I3.7.3P but also to the analysis method I3.7.4p . 

• Flag 2 is not used. We didn’t go so far as to rerun the full detection and localization 
procedure (§ I3.ip with the alternative diffuse model. Assessing the changes in source 
characteristics is normally enough. 

• For Flag 4, we lowered the threshold for flagging the source-to-background ratio to 10%, 
recognizing that the uncertainties in the interstellar emission model are now reduced 
(App. 1^. 

• We reinstated Flag 7 (comparison between pointlike and gtfindsrc localizations) which 
was not used in 2FGL because of an inconsistency in the unbinned likelihood results. 
It indicates sources for which the source locations derived from pointlike (§ I3.1.3P 
and gtfindsrc are inconsistent at the 95% conhdence level, gtfindsrc was applied 
only above 3 GeV due to computing time constraints. This is appropriate for most 
sources (because the PSF is much better at high energy) but does not allow testing 
the localization of soft sources. 

• Flag 8 has been merged into Flag 9. Both tested localization reliability. 

• Flag 11 is deprecated because we put in place an explicit time-dependent model for 
the Sun and Moon emission (f 12.3p . 


4. The 3FGL Catalog 

We present a basic description of the 3FGL catalog in § 14.11 including a listing of the 
main table contents and some of the primary properties of the sources in the catalog. We 
present a detailed comparison of the 3FGL catalog with the 2FGL catalog in § 14.21 


4.1. Catalog Description 

Table H] is the catalog, with information for each of the 3033 source^; see Table |5] for 
descriptions of the columns. The source designation is 3FGL JHHMM.m+DDMM where the 3 
indicates that this is the third LAT catalog, FGL represents Fermi Gamma-ray LAT. Sources 
close to the Galactic ridge and some nearby interstellar cloud complexes are assigned names 


^^Table|4]has 3034 entries because the PWN component of the Crab nebula is represented by two cospatial 
sources G 13.311 . 
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Fig. 15.— Full sky map (top) and blow-up of the inner Galactic region (bottom) showing 
sources by source class (see Table [6]). All AGN classes are plotted with the same symbol for 
simplicity. 
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of the form 3FGL JHHMM.m+DDMMc, where the c indicates that caution should be used in 
interpreting or analyzing these sources. Errors in the model of interstellar diffuse emission, 
or an unusually high density of sources, are likely to affect the measured properties or even 
existence of these 78 sources (see I3.8p . In addition a set of analysis flags has been dehned 
to indicate sources with unusual or potentially problematic characteristics (see 13.91) . The 
c designator is encoded as one of these flags. An additional 572 sources have one or more of 
the other analysis flags set. The 25 sources that were modeled as extended for 3FGL 13.41) 
are singled out by an e appended to their names. 

The designations of the classes that we use to categorize the 3FGL sources are listed 
in Table [ 6 ] along with the numbers of sources assigned to each class. Figure [15] illustrates 
where the source classes are in the sky. We distinguish between associated and identihed 
sources, with associations depending primarily on close positional correspondence (see § 15.2p 
and identihcations requiring measurement of correlated variability at other wavelengths or 
characterization of the 3FGL source by its angular extent (see 11 15.ip . In the cases of multiple 
associations with a 3FGL source, we adopt the single association that is statistically most 
likely to be true if it is above the conhdence threshold (see § 15.2p . Sources associated 
with SNRs are often also associated with PWNe and pulsars, and the SNRs themselves are 
often not point-like. We do not attempt to distinguish among the possible classihcations 
and instead list in Table [7| plausible associations of each class for unidentihed 3FGL sources 
found to be positionally associated with SNrJ^. The Crab pulsar and PWN are represented 
by a total of three entries, two of which (designated i and s) represent spectral components 
of the PWN (see 11 15.ip . We consider these three entries to represent two sources. 

The photon flux for 1-100 GeV (T 35 ) and the energy flux for 100 MeV to 100 GeV 
in Table 01 are evaluated from the £t to the full band (see § 13.5p . We do not present the 
integrated photon flux for 100 MeV to 100 GeV (see § 13.5p . Table [S] presents the fluxes in 
individual bands as defined in § 13.51 


^"^Four sources positionally associated with SNRs were also found to be associated with blazars. We can¬ 
not quantitatively compare association probabilities between the blazar and the (spatially extended) SNR 
classes. In the 3FGL catalog, we list only the blazar associations for them. The sources and SNR asso¬ 
ciations are 3FGL J0217.3-b6209 (G137.2-h01.3), 3FGL J0223.5-b6313 (G132.7-b01.3), 3FGL J0526.0-b4253 
(G166.0-1-04.3), and 3FGL J0215.6-b3709 (G074.9-b01.2). 








Table 4. LAT 4-year Catalog 


Name 3FGL 

R.A. 

Decl. 

1 

b 

di 

02 </> 

CT 

J^35 

AF 35 

•S '25 

AS 25 

r25 

Ar25 

Mod 

Var 

Flags 

7 -ray Assoc. TeV 

Glass® 

ID or Assoc. 

JOOOO. 1+6545 

0.038 

65.752 

117.694 

3.403 

0.102 

0.078 41 

6.8 

1.0 

0.2 

13.6 

2.1 

2.41 

0.08 

PL 


3 

2FGL J2359.6+6543C ■ ■ ■ 



JOOOO.2-3738 

0.061 

-37.648 

345.411 

-74.947 

0.073 

0.068 -89 

5.1 

0.2 

0.1 

2.4 

0.7 

1.87 

0.18 

PL 






JOOOl.0+6314 

0.254 

63.244 

117.293 

0.926 

0.248 

0.160 -65 

6.2 

0.6 

0.1 

13.0 

1.9 

2.73 

0.11 

PL 


3,4,5 

2FGL J2358.9+6325 ■ ■ ■ 

spp 


JOOOl.2-0748 

0.321 

-7.816 

89.022 

-67.324 

0.082 

0.070 -19 

11.3 

0.7 

0.1 

7.8 

0.9 

2.15 

0.09 

PL 



2FGL JOOOO.9-0748 ■ ■ ■ 

IFGL J0000.9-0745 

bll 

PMN JOOOl-0746 

JOOOl.4+2120 

0.361 

21.338 

107.665 

-40.047 

0.211 

0.188 -33 

11.4 

0.3 

0.1 

8.1 

0.8 

2.78 


LP 

T 


3EG J2359+2041 

fsrq 

TXS 2358+209 

JOOOl.6+3535 

0.404 

35.590 

111.661 

-26.188 

0.213 

0.167 8 

4.2 

0.3 

0.1 

3.4 

0.8 

2.35 

0.19 

PL 


4 




J0002.0-6722 

0.524 

-67.370 

310.139 

-49.062 

0.102 

0.086 69 

5.9 

0.3 

0.1 

3.3 

0.8 

1.95 

0.16 

PL 






J0002.2-4152 

0.562 

-41.883 

334.070 

-72.143 

0.217 

0.140 68 

5.2 

0.3 

0.1 

3.0 

0.7 

2.09 

0.19 

PL 



2FGL JOOOl.7-4159 ■ ■ ■ 

bcu 

IRXS 

J000135.5-415519 


















IFGL JOOOl.9-4158 



J0002.6+6218 

0.674 

62.301 

117.302 

-0.037 

0.061 

0.054 -55 

18.0 

2.8 

0.2 

18.4 

1.7 

2.35 


LP 



2FGL J0002.7+6220 ■ ■ ■ 



J0003.2-5246 

0.815 

-52.777 

318.976 

-62.825 

0.070 

0.061 -44 

5.7 

0.3 

0.1 

3.0 

0.8 

1.90 

0.17 

PL 




bcu 

RBS 0006 

J0003.4+3100 

0.858 

31.008 

110.964 

-30.745 

0.181 

0.163 13 

6.3 

0.3 

0.1 

4.9 

0.8 

2.55 

0.13 

PL 






J0003.5+5721 

0.890 

57.360 

116.486 

-4.912 

0.089 

0.072 1 

5.4 

0.5 

0.1 

5.4 

1.1 

2.18 

0.13 

PL 





1 

-- 


^See Table [6| for class designators. I 


Note. 

content. 


This table is published in its entirety in the electronic edition of the Astrophysical Journal Supplements. A portion is shown here for guidance regarding its form and 
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Table 5. LAT Third Catalog Description 


Column 


Description 


Name 


R.A. 

Decl. 

I 

b 

0i 

O 2 

<P 

<j 

^35 

AFss 

5'25 

AS'25 

r25 

Ar25 

Mod. 

Var. 

Flags 

7-ray Assoc. 

TeV 

Class 

ID or Assoc. 


3FGL JHHMM.m+DDMM[c/e/i/s], constructed according to lAU Specifications for Nomenclature; m is decimal 

minutes of R.A.; in the name, R.A. and Decl. are truncated at 0.1 decimal minutes and 1', respectively; 

c indicates that based on the region of the sky the source is considered to be potentially confused 

with Galactic diffuse emission; e indicates a source that was modeled as spatially extended (see ^ 13.41 1: 

the two spectral components of the Crab PWN are designated i and s 

Right Ascension, J2000, deg, 3 decimal places 

Declination, J2000, deg, 3 decimal places 

Galactic Longitude, deg, 3 decimal places 

Galactic Latitude, deg, 3 decimal places 

Semimajor radius of 95% confidence region, deg, 3 decimal places 

Semiminor radius of 95% confidence region, deg, 3 decimal places 

Position angle of 95% confidence region, deg. East of North, 0 decimal places 

Significance derived from likelihood Test Statistic for 100 MeV-300 GeV analysis, 1 decimal place 

Photon flux for 1 GeV—100 GeV, 10~^ ph cm“^ summed over 3 bands, 1 decimal place 

1(7 uncertainty on F 35 , same units and precision 

Energy flux for 100 MeV—100 GeV, 10“^^ erg cm“^ s~^, from power-law fit, 1 decimal place 
1(7 uncertainty on 525 , same units and precision 

Photon number power-law index, 100 MeV-100 GeV, 2 decimal places 

1(7 uncertainty of photon number power-law index, 100 MeV-100 GeV, 2 decimal places 

PL indicates power-law fit to the energy spectrum; LP indicates log-parabola fit to the energy spectrum; 

EC indicates power-law with exponential cutoff fit to the energy spectrum 
T indicates < 1% chance of being a steady source; see note in text 
See Table [ 3 ] for definitions of the flag numbers 

Positional associations with OFGL, IFGL, 2FGL, 3EG, EGR, or lAGL sources 

Positional association with a TeVGat source, P for unresolved angular size, E for extended 

Like ‘ID’ in 3EG catalog, but with more detail (see Table[ 6 ll. Gapital letters indicate firm identifications; 

lower-case letters indicate associations 

Designator of identified or associated source 
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Table 6. LAT 3FGL Source Classes 


Description 

Identified 


Associated 


Designator Number 

Designator 

Number 

Pulsar, identified by pulsations 

Pulsar, no pulsations seen in LAT yet 

PSR 

143 

psr 

24 

Pulsar wind nebula 

PWN 

9 

pwn 

2 

Supernova remnant 

SNR 

12 

snr 

11 

Supernova remnant / Pulsar wind nebula 



spp 

49 

Globular cluster 

GLC 

0 

glc 

15 

High-mass binary 

HMB 

3 

hmb 

0 

Binary 

BIN 

1 

bin 

0 

Nova 

NOV 

1 

nov 

0 

Star-forming region 

SFR 

1 

sir 

0 

Compact Steep Spectrum Quasar 

CSS 

0 

CSS 

1 

BL Lac type of blazar 

BLL 

18 

bll 

642 

FSRQ type of blazar 

FSRQ 

38 

fsrq 

446 

Non-blazar active galaxy 

AGN 

0 

agn 

3 

Radio galaxy 

RDG 

3 

rdg 

12 

Seyfert galaxy 

SEY 

0 

sey 

1 

Blazar candidate of uncertain type 

BCU 

5 

bcu 

568 

Normal galaxy (or part) 

GAL 

2 

gal 

1 

Starburst galaxy 

SBG 

0 

sbg 

4 

Narrow line Seyfert 1 

NLSYl 

2 

nlsyl 

3 

Soft spectrum radio quasar 

SSRQ 

0 

ssrq 

3 

Total 


238 


1785 

Unassociated 




1010 


Note. — The designation ‘spp’ indicates potential association with SNR or PWN (see 
Table O . Designations shown in capital letters are firm identifications; lower case letters 
indicate associations. In the case of AGN, many of the associations have high confidence. 
Among the pulsars, those with names beginning with LAT were discovered with the LAT. 







Table 7. Potential Associations for Sources Near SNRs 


Name 3FGL 

SNR name 

PWN name 

TeV name 

Common name 

JOOOl.0+6314 

G116.5+01.1 




J0128.4+6257 

G127.1+00.5 



R5 

J0220.1+6202C 

G132.7+01.3 



HB3 

J0224.0+6235 

G132.7+01.3 



HB3 

J0500.3+5237 

G156.2+05.7 




J0610.6+1728 

G192.8-01.1 



PKS 0607+17 

J0631.6+0644 

G205.5+00.5 



Monoceros Loop 

J0640.9+0752 

G205.5+00.5 



Monoceros Loop 

30838.1-4615=' 

G263.9-03.3 



Vela 

J0839.1-4739 

G263.9-03.3 



Vela 

J0843.1-4546 

G263.9-03.3 



Vela 

JllOl.9-6053 

G290.1-00.8 



MSH 11-61A 

Jllll.9-6038 

G291.0-00.1 

G291.0-0.1 


MSH 11-62 

J1209.1-5224 

G296.5+10.0 



PKS 1209-51 

J1212.2-6251 

G298.5-00.3 




J1214.0-6236 

G298.6-00.0 




J1345.1-6224 

G308.8-00.1 




J1441.5-5955C 

G316.3-00.0 



MSH 14-57 

J1549.1-5347C 

G327.4+00.4 




J1551.1-5610 

G326.3-01.8 



Kes 25 

J1552.9-5610 

J1615.3-5146e 

G326.3-01.8 

HESS 31614-518 


Kes 25 

J1628.9-4852 

G335.2+00.1 




J1636.2-4709C 

G337.2+00.1 


HESS 31634-472 


31638.6-4654= 

G337.8-00.1 



Kes 41 

31640.4-4634C 

G338.3-00.0 


HESS 31640-465 


31641.1-4619C 

G338.5+00.1 


HESS 31641-463 


31645.9-5420 

G332.5-05.6 




31722.9-4529 

G343.0-06.0 



RGW 114 

31725.1-2832 

G358.0+03.8 




31728.0-4606 

G343.0-06.0 



RGW 114 

31729.5-2824 

G358.0+03.8 




31737.3-3214C 

G356.3-00.3 




31745.1-3011 

G359.1-00.5 


HESS 31745-303 


31745.6-2859C 

GOOO.0+00.0 

G359.95-0.04 

Galactic Centre 

Sgr A East 

31810.1-1910 

GOll.1+00.1 


HESS 31809-193 


31811.3-1927C 

GOll.2-00.3 


HESS 31809-193 


31817.2-1739 

G013.3-01.3 




31818.7-1528= 

G015.4+00.1 




31828.4-1121= 

G020.0-00.2 




31829.7-1304 

G018.9-01.1 




31833.9-0711= 

G024.7+00.6 




31834.6-0659 

G024.7+00.6 




31840.1-0412 

G027.8+00.6 




31915.9+1112 

G045.7-00.4 




31951.6+2926 

G065.7+01.2 
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4.2. Comparison with OFGL, IFGL, 2FGL and IFHL 

4-2.1. General comparison 

We compare the 3FGL with previous catalogs released by the LAT collaboration. These 
are listed in Tabled 

Associations among 3FGL and 0FGL/1FGL/2FGL and IFHL sources are based on the 
following relation: 

A < dx,a = + ^x,3FGL 

where A is the angular distance between the sources, is dehned in terms of the semi-major 
axis of the x% conhdence error ellipse for the position of each source, e.g., the 95% conhdence 
error for the automatic source association procedure (§ 15.21) and ‘a’ is alternatively OFGL, 
IFGL, 2FGL, and IFHL. In total, 1720 3FGL sources were automatically associated with 
entries in either the OFGL, IFGL, 2FGL, or IFHL catalogs. The statistics of the association 
results is reported in Table [131 

In the 3FGL analysis the spectral hts are made using power-law, power-law with an 
exponential cutoff, or log-parabola models 13.3p . For each 2FGL and 3FGL source we also 
evaluated the spectral index (T) of the best power-law £t and this enables a comparison of 
the spectral characteristics of the IFGL, 2FGL, and 3FGL sources. Figure HE] shows the 
distributions of the power-law indices of the sources at high Galactic latitude and only those 
with a power-law spectral type in the IFGL, 2FGL and 3FGL catalogs, to avoid possible 
confusion from more complex features. The three distributions are very similar, with an 
average Tifgl = 2.23 ± 0.01, average r 2 FGL = 2.21 ± 0.01, and average Tsfgl = 2.19 ± 0.01. 
However, the peaks of the three distributions are not exactly coincident; also, they have 


Table 7—Gontinued 


Name 3FGL 

SNR name PWN name 

TeV name Common name 

J2014.4+3606 

G073.9+00.9 ■ ■ ■ 


J2022.2+3840 

G076.9+01.0 ■ ■ ■ 


J2225.8+6045 

G106.3+02.7 ■■■ 

G106.3+2.7 


^These sources have been found to be significantly variable, i.e., Variability.Index > 72.44 I3.6|l . which would be unex¬ 

pected for physical associations with SNRs or PWNe. 


Note. — These sources are classified as spp in Table [4| They may be pulsars rather than the SNR or PWN named. Four 
additional 3FGL sources are associated with both an SNR and a blazar. For these the catalog lists the blazar associations; see 
text. 
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Table 8. LAT 4-year Catalog: Spectral Information 



0.1-0.3 GeV 

0.3-1 GeV 

1-3 GeV 

3-10 GeV 


10-100 GeV 

Name 3FGL 


v/rsT 

F2^ 

v/rsj 


VT^ 

Fi’^ V 

T54 

Fs^ 

VTSi 

JOOOO. 1+6545 

J--o1_o.84 

2.2 

u.oy_Q 14 

5.3 

1 24+0-^^ 
-^•^^-0.23 

6.0 

0-58l°;« 

1.5 

0.28t°;« 

2.4 

JOOOO.2-3738 


0.0 

U-Ui_Q Ql 

0.6 

o-i54g;g^ 

2.9 

0.75t°oil 

4.2 

o-ietgil 

3.5 

JOOOl.0+6314 

9 Ql +0-74 

4.0 

-0.11 

4.3 

0.30l°;« 

1.7 

1 lfi+0-48 
J--10_o.42 

3.3 

0 02+0-^^ 
C.UZ-o 02 

0.4 

JOOOl.2-0748 

0 41 

1.6 

0 21 

U-^l_0.04 

5.3 

0 52+0’^^ 

7.0 

9 n-b0.5i 

^•4J-_o.45 

8.0 

n 10+0-15 

'J-1C_0.08 

1.6 

JOOOl.4+2120 

. c-9+0.24 

i.OZ_Q 24 

6.8 

U.OD_0 05 

7.6 

0 ‘:{f5+0-10 

C.OD_o.o9 

4.5 

0 oo+o-^o 
U.UU_o.oo 

0.0 

0 nn+0-^^ 
C.UU-o.oo 

0.0 

JOOOl.6+3535 

0.95lg;2? 

3.6 

U-'J^_0.04 

0.9 

Q ly-l-O.OO 

'J-lf _0.08 

2.5 

0.40l°Ji 

2.0 

0.23l°;« 

2.3 

J0002.0-6722 

0-054g:g| 

0.2 

U.UD_o 03 

2.0 

0 20+0 08 
C.^U_o.o7 

3.4 

0.83l°i2 

4.8 

0.28l°;« 

3.1 

J0002.2-4152 

U.OO_Q 22 

2.5 

0.06l“;°| 

2.1 

0 09+0-07 

u.uy_o.o6 

1.5 

0.98l°;|® 

4.7 

0.26l°;« 

2.6 

J0002.6+6218 

i.UO_Q 4Q 

2.6 

0.80t°;°® 

9.4 

2.44t°i| 

13.3 

3 5«+0-69 

'J.OO_o.64 

7.4 

o.ool°;J® 

0.0 

J0003.2-5246 

0.66lg;|» 

1.8 

U.UU_o.oo 

0.0 

0.l8l°;“« 

3.1 

0 70+0-36 
' ^-0.28 

3.7 

0.30l°;« 

4.5 

J0003.4+3100 

0.86t“:=^ 

3.4 

0 15+004 

'J-J-O_o.o4 

3.6 

o.364°;J“ 

4.9 

0 20+0-^0 
C.^U_o.i9 

1.1 

0 01 +0.14 

C.Ui_0.0l 

0.1 

J0003.5+5721 

U.UU_Q QQ 

0.0 

0 23+0 07 

u.zo_o.o7 

3.5 

o-494g;ll 

4.0 

0 02+0-41 
'J-y^_0.35 

3.3 

0.2itVlt 

3.0 

J0003.8-1151 

U.Ui_o.oi 

0.0 

0.054°:°^ 

1.5 

n 9 q-|- 0.08 

^•'^‘^-0.07 

3.8 

0 47+0-30 
-0.24 

2.5 

0.l8l°;l® 

2.5 

J0004.2+6757 

U.Oy_Q 59 

0.5 

0.56l“;« 

4.4 

0 59+0-10 

*J-oy_o.i8 

3.4 

1 15+0-47 
J---LO-0.42 

3.2 

0.09l°;“ 

1.4 

J0004.2+0843 

f\ r)q + 0.25 

U.zo_o.23 

0.8 

0.03t“:°4 

0.7 

0 12+009 

LI-J-^_0.08 

1.5 

1 lfi+0-41 
4-19_o.35 

5.3 

0-08tg;Jg 

2.0 


Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal Supplements. A portion 
is shown here for guidance regarding its form and content. 

^In units of 10“^ photons cm“^ s“^ 

'^In units of 10~^ photons cm“^ s~^ 

^In units of 10“^^ photons cm~^ s“^ 


Table 9. Statistics of Sources in LAT Catalogs 


Gategory 

OFGL’* 

IFGL 

2FGL 

IFHL^ 

3FGL 

Total 

205 

1451 

1873 

514 

3033 

High-Latitude sources 

132 

1043 

1319 

399 

2193 

Low-Latitude sources 

73 

408 

554 

115 

841 

‘Lost’ sources('’’bn 3FGL 

12 

310 

300 

17 

- 


^OFGL, the LAT Bright Source List, has a lower energy limit of 200 MeV 
and a significance threshold TS > 100. 

^IFHL is a catalog for the energy range >10 GeV. 

^Sources without a counterpart in 3FGL catalog, at the level of overlap¬ 
ping 95% source location confidence contours. These sources are discussed in 
Table 1101 
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Fig. 16.— Distributions of the spectral index for the high-latitude sources (|6| > 10°) in 
IFGL (1043 sources, blue line), 2FGL (1173 sources, red line) 3FGL (1960 sources, black 
line) catalogs. 2FGL and 3FGL samples include only power-law spectrum type. 

different skewnesses. The small differences in the power-law index distributions could be 
related to slightly different systematic uncertainties in the effective area between the in¬ 
strument response functions P7REP_SOURGE_V15, P7SOURGE_V6 and P6_V3_DIFFUSE 
used respectively for 3FGL, 2FGL, and IFGL. 

We have compared the distribution of the 95% conhdence error radii of the IFGL, 
2FGL, and 3FGL sources at high Galactic latitude. The distribution of 95% conhdence 
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95% confidence error radius (arcmin) 


Fig. 17.— Distributions of the 95% confidence error radii for high-latitude sources (|6| > 10°) 
with 25 < TS < 100 in IFGL (blue line), 2FGL (red line) and 3FGL (black line), illustrating 
the improvement of localizations for sources of equivalent detection signihcances. 

error radius for those sources with 25 < TS < 100 in any of the IFGL, 2FGL, and 3FGL 
catalogs (Figure [T7]) shows the localization improvement for a given range of source detection 
signihcances. We evaluated the 95% conhdence error radius as the geometric mean of the 
semi-major and semi-minor axes of the 95% conhdence error ellipse. 

Figure [TS] shows the energy hux distribution in IFGL, 2FGL, and 3FGL. Comparing 
the current hux threshold with those published in previous LAT Catalog papers we see that 
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Energy Flux 100 MeV — 100 GeV (erg/cm^/s) 


Fig. 18.— Distributions of the energy flux for IFGL (blue line), 2FGL (red line), and 3FGL 
(black line) sources at high Galactic latitude {\b\ > 10°). 

in 3FGL the threshold is down to ~ 3 x 10“^^ erg cm“^ s“^, from ~ 5 x 10“^^ erg cm“^ 
s“^ in 2FGL and ~ 8 x 10“^^ erg cm“^ s“^ in IFGL. Above that flux the 2FGL and 3FGL 
distributions are entirely compatible. 

However the IFGL distribution shows a distinct bump between 1 and 2 x 10“^^ erg cm“^ 
s“^. That accumulation of fluxes was clearly incorrect. We attribute it primarily to overes¬ 
timating signihcances and fluxes due to the unbinned likelihood bias in the IFGL analysis, 
and also to the less accurate procedure then used to extract source flux (see discussion in 
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the 2FGL paper). 


4 . 2 . 2 . Comparison of individual sources 



Fig. 19.— Distribution of the differences Fsfgl — Tifgl (blue line) and Fsfgl — r 2 FGL (red 
line) for the 621 sources at high latitude (|6| > 10°) in common among the IFGL, 2FGL 
and 3FGL catalogs. For the 2FGL and 3FGL samples only power-law spectrum type sources 
have been considered. 

Figure [19] shows the distribution of the differences Fsfgl — r 2 FGL and Fsfgl — Tifgl 
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for the 621 high-latitude sources with power-law spectrum type in common among the three 
catalogs. The average of the 3FGL - 2FGL distribution is 0.04±0.01, with the 3FGL sources 
slightly softer than the 2FGL ones, while the average of the 3FGL - IFGL distribution is 
—0.04 ± 0.01, with the 3FGL sources slightly harder than the IFGL ones. 



95% confidence error radius (arcmin) 


Fig. 20.— Distributions of the 95% confidence error radius for high-latitude sources (|6| > 
10°) in common among IFGL (blue line), 2FGL (red line), and 3FGL (black line), illustrating 
the improvement of localizations for sources of equivalent detection significances. 

When comparing the distribution of 95% conhdence error radius for the sources in 
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common among all the LAT catalogs, we see that for 3FGL this parameter extends to lower 
values than for the earlier catalogs, showing that the localization has improved, thanks to 
improvements in the 3FGL analysis and increased statistics over the longer integration period 
for 3FGL (Figure 1201). 


4 . 2 . 3 . Possible causes for losing sources 

In the remainder of this section we describe a variety of reasons why the ‘lost’ OFGL, 
IFGL, 2FGL, and IFHL sources might not appear in the 3FGL catalog. Table [TUI shows the 
statistics of the ‘lost’ sources. 

We have also produced tables with all the ‘lost’ sources for each previous LAT catalog. 
The first rows of the ‘lost’ 7 -ray source table for the 2 FGL catalog are listed in Table [TTl 
only reported for guidance^ In the last columns of Table [TTl we assigned to each source one 
or more flags corresponding to possible causes for it to be lost and which we will discuss in 
the following paragraphs. In many cases, no one reason can be singled out. 

Lost sources from previous LAT catalogs are in general equally distributed over all 
latitudes, with a slight excess at low latitudes for 2FGL ‘lost’ sources. In fact about 10% 
of the 2FGL ‘lost’ sources are at low Galactic latitude compared to a 6 % of high-latitude 
‘lost’ 2FGL sources. We remind the reader that at low latitudes the Galactic diffuse emission 
is most intense and improvements in the model for the diffuse emission would be expected 
to have the most influence (§ 12.3p . The sources in common among 3FGL and the previous 
LAT catalogs are primarily outside the Galactic plane, as are the sources newly detected 
in 3FGL. Most of the ‘lost’ sources were also listed as unassociated in the previous FGL 
catalogs. Among the former associated ‘lost’ sources, most of them were associated with 
AGN and a few with pulsars. For sources of the AGN type their absence from the 3FGL 
catalog can be due to their intrinsic variability. A faint source which flared during the hrst 
year, allowing it to be detected in OFGL, can be diluted and become undetectable in a longer 
time interval. 

Most of the ‘lost’ sources have analysis flags or the c designator in IFGL and 2FGL 
names, indicating that these sources were already flagged as influenced by the diffuse emission 
and recognized as potentially problematic or possibly spurious. 

Some other IFGL, 2FGL, and IFHL sources do not have counterparts in the 3FGL 


^^The full table of lost 2FGL sources and similar tables for lost OFGL, lost IFGL, and lost IFHL sources 
are available only in the electronic version. 




Table 10. Statistics of ‘Lost’ OFGL, IFGL, 2FGL, and IFHL Sources 



OFGL not in 3FGL 

IFGL not in 3FGL 

2FGL not in 3FGL 

IFHL not in 3FGL 

All 

12 

310 

300 

17 

With flags 

- 

131 

211 

- 

Name-FGL c 

- 

104 

87 

- 

AGN 

1 

22 

27 

1 

PSR 

0 

1 

3 

0 

Unassociated 

11 

264 

234 

16 

Within 1° of a 3FGL e 

3 

27 

33 

4 

sources in other FGL catalogs 

OFGL 

- 

5 

5 

0 

IFGL 

4 

- 

56 

1 

2FGL 

3 

67 

- 

1 

IFHL 

0 

2 

8 

- 

Not in any other Fermi catalog 

7 

237 

237 

15 


indicates that based on the region of the sky the source is considered to be potentially confused with Galactic diffuse emission, 
indicates a source that was modeled as spatially extended. 
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catalog because they have been resolved into two or more 3FGL sources or candidate source 
seeds. We flag them with ‘S’ (split) in the ‘Flag’ column of Table fTTl In some cases only one 
of the seeds reached TS > 25 and so was included in the 3FGL list. 

Several other possible causes of ‘lost’ sources are evident: (1) the 3FGL y-ray centroid 
has shifted with respect to the previous FGL catalogs, preventing the matching; (2) statisti¬ 
cal threshold effects, i.e. their TS has dropped below 25. Additional considerations include 
variability and (generally small) effects from the different event selections used for the anal¬ 
yses (P7REP_SOURGE_V15 for 3FGL, P7GLEAN_V6 for IFHL, P7SOURGE_V6 for 2FGL 
and P6_V3_DIFFUSE for OFGL); different Galactic diffuse emission models; different analy¬ 
sis procedures (unbinned likelihood analysis for OFGL and IFGL, binned likelihood analysis 
for 2FGL and IFHL, and a combination of binned and unbinned for 3FGL). We analyze 
those causes in more detail for 2FGL in § 14.2.41 We stress that these differences are often 
not negligible. 

A comparison of the source significances of the ‘lost’ sources with those in the 3FGL 
catalog shows that (Figure 1211) in the latter we have not lost highly-signiflcant sources. The 
peaks of the source significance distributions for all the sources of the FGL catalogs (not 
shown in the Figure [2T]) have shifted from 4-6cr for IFGL to 4-5(j for 2FGL and 3FGL. 

The power-law indices of high-Galactic latitude (|6| > 10°) ‘lost’ sources with power-law 
spectral type tend to be softer than average for their catalogs (Figure [22il . 

The numbers of associated sources among the OFGL, IFGL, and 2FGL catalogs, and 
the 3FGL catalog do depend on the criteria used to define spatial coincidence (Eq. 0]). The 
numbers of IFGL - 3FGL, 2FGL - 3FGL and IFHL - 3FGL associated sources increase if we 
use A < dgg.g as association criteriofl The 193 additional associations (listed in Table M 
and corresponding OFGL, IFGL, IFHL tables in the column ‘3FGL (A < dgg.g)’) represent 
about 5% of the OFGL, IFGL, 2FGL, and IFHL sources, as expected when passing from 
dgs to dgg.g. Furthermore, the improved model of the Galactic diffuse emission 12.3p used 
to build the 3FGL catalog together with the expected increase of the signal-to-noise ratio 
due to the use of 48 months of data, allowed us to obtain better localizations of the sources 
at positions that might be outside the 95% confidence error regions reported in the OFGL, 
IFGL, or 2FGL catalogs. Indeed, about half of the 193 additional associations concern 
sources located along the Galactic plane. Also, in the IFGL catalog the positions of sources 
associated with the LAT-detected pulsars and X-ray binaries are the high-precision positions 
of the identified counterparts. (These sources can be easily recognized because they have 


Assuming a Rayleigh distribution for the source angular separations, dgg.g is evaluated using 
%9.9 = 1-52 0g5. 
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Source significance (cr) 


Fig. 21.— Distributions of the significances of ‘lost’ IFGL and ‘lost’ 2FGL sources compared 
to the IFGL and 2FGL sources which are associated to 3FGL sources. All sources at high 
Galactic latitudes (|6| > 10°) are included. (2FGL sources associated to 3FGL sources: solid 
red line, ‘lost’ 2FGL sources: dashed red line, IFGL sources associated to 3FGL sources: 
solid blue line, ‘lost’ IFGL sources: dashed blue line). 

null values in the localization parameters reported in the IFGL catalog.) Not all of these 
associations appear in the 3FGL catalog because they cannot be associated using dg^, but 
those that can be associated using dgg,g are listed in Table [TT] (and corresponding OFGL, 
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Power-law index 


Fig. 22.— Distributions of the power-law index of the IFGL sources (blue solid line) and 
2FGL sources (red solid line) in the 3FGL catalog and of the ‘lost’ IFGL sources (blue 
dashed line) and 2FGL sources (red dashed line). All samples include only high-latitude 
sources (|6| > 10'’). 

IFGL, and IFHL tables). 

To study a possible reason for OFGL, IFGL, 2FGL, and IFHL sources to disappear in 
the 3FGL catalog, we have compared the TS they had when published in their respective 
catalogs with their values in the 3FGL pointlike analysis. The 3FGL catalog was built, in 
fact, starting from 4003 seeds with TS* > 10 in the pointlike analysis (§ 1^. The hnal gtlike 
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analysis, which did not change the positions of the seeds, resnlted in the 3033 sources with 
TS > 25 that make up the 3FGL catalog. Therefore, possibly many seeds did not reach 
the threshold but can be associated with OFGL, IFGL, 2FGL, and IFHL sources (using 
A < dgg.g). These sources, marked with ‘T’ (for True’) in the ‘3FGL Seed’ column of the 
Table mi can be considered to be previously confirmed sources whose significance dropped 
below the threshold, either as a result of time variability, changes in the model or in the 
catalog analysis procedure for Galactic diffuse emission. Finally we looked for those ‘lost’ 
sources whose distances from an extended 3FGL source are less than 1°, and these are flagged 
with ‘E’ in the column ‘Flag’ of Table fTTl 


4.2.4- Step-by-step from 2FGL to 3FGL 

In order to understand the improvements of the 3FGL analysis with respect to 2FGL, 
we have considered the effects of changing the analysis, the data set, and the diffuse emission 
model without changing the time range (i.e., leaving it as two years). To that end we started 
with the 2FGL catalog and changed each of those three elements in sequence and compared 
each intermediate result with the previous one. 

• The main difference between the analyses is the Front/Back handling (§ 13.2p . The 
comparison showed that using identical isotropic diffuse spectra for Front and Back 
events in 2FGL resulted in underestimating the low-energy fluxes of high-latitude 
sources. As a consequence, the corrected analysis leads to larger TS values, higher 
photon fluxes, softer spectra, and smaller curvatures than in 2FGL. The effects are 
small on the scale of individual sources but collectively obvious. Quantitatively, the 
average difference in spectral index induced by this change was measured to be -1-0.05. 
Because that effect is due to the background, it is at the same level in a units (~ O.da) 
for faint and bright sources. 

• Ghanging data from Pass 7 (2FGL) to Pass 7 reprocessed (3FGL) results in some¬ 
what larger TS", harder sources and more curved spectra (but no change of integral 
flux on average). The average difference in spectral index is —0.03. This goes in the 
opposite direction to (and therefore partly offsets) the difference due to the separate 
Front/Back handling. However the dependence on flux is not the same. The repro¬ 
cessing affects essentially all spectral indices and curvatures equally in absolute terms. 

• Finally, changing the model for Galactic diffuse emission from gal_2yearp7v6_v0 used 
in 2FGL to gll_iem_v05_revl results in smaller TF, lower fluxes and less curved 
spectra (but no change of spectral index on average). Like the first point above, this 
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background-related effect is smaller in absolute (curvature) or relative (flux) terms for 
brighter sources. 

In conclusion, to first order the resulting net changes are not very large, consistent with the 
general comparison between 3FGL and 2FGL at the beginning of this section. The 3FGL 
sources tend to be less curved than 2FGL ones. In particular, there are fewer pathological 
very strongly curved sources (with /? = 1 and Flag 12 set) in 3FGL (41) than in 2FGL (64) 
even though there are more LogParabola spectra in 3FGL (395) than in 2FGL (336) because 
of the better statistics. 


5. Source Association and Identification 
5.1. Firm Identifications 

As with the 2FGL and earlier LAT catalogs, we retain the distinction between associa¬ 
tions and firm identihcations. Although many associations that we list between LAT sources 
and potential counterparts at other wavelengths, particularly those for AGN, have very high 
probability of being true, a firm identification, shown in the catalog by capitals in the Glass 
column in Table O is based on one of three criteria: 


1. Periodic Variability. Pulsars are the largest class in this category. All PSR labels 
indicate that pulsed 7 rays have been seen from the source with a probability of the 
periodicity occurring by chance of less than 10“®. Pulsars detected in blind searches of 
LAT data are indicated as ‘LAT PSR’ in the ‘ID or Assoc.’ column of Table H the other 
PSR detections are based on folding with radio or X-ray ephemerides (see lAbdo et ah 


20131 ). A similar chance probability requirement applies to the other set of periodic 


sources, t he high-mass binar ies (HMB). Three of these are i ncluded in the catalog: LSI 
+61 303 flAbdo et aLl[2nn^ . LS 5039 flAbdo et al.ll2009eh . and IFGL J1018.6-5856 


(iGorbet et al.ll201ll) . Although not quite meeting the same cha nce probability, another 


binar y (BIN) is included as an identification: Eta Garinae (IReitberger et al.l 12012 


20141 1. 


2 . Spatial Morphology. Spatially extended sources whose morphology can be related to 
extent seen at other wavelengths include SNR, PWNe, and galaxies, as described in 
§ 13.41 The Gentaurus A lobes and core are both marked as identified, because they 
are part of the same extended source, although the core itself does not show spatial 
extent. Although individual molecular clouds could in principle be included in this 





























Table 11. 2FGL Sources Not in the 3FGL Gatalog 


2FGL 

2FGL Assoc. 

1 ( 0 .) 

(deg) 

(deg) 

(deg) 

o-(“) 

Spec/“^ 

Index 

Var(“) 
p > 90% 

3FGL('’) 

(A < 0(99.9) 

SFGL^'^) 

( 0 ( 99.9 < A < 1°) 

(deg) 

A /dgg ,9 

3FGL(®) FlagsW) 
Seed 

J0004.2+2208 


108.732 

-39.430 

0.194 

5.4 

2.49 

F 





T 


JOOll.3+0054 

PMN JOOll+0058 

102.317 

-60.352 

0.223 

6.2 

2.43 

T 






F 

J0013.8+1907 

GB6 J0013+1910 

110.786 

-42.858 

0.160 

4.1 

2.06 

T 







J0014.3-0509 


99.374 

-66.312 

0.169 

4.2 

2.45 

F 

J0014.3-0455 


0.322 

0.725 

T 

F 

J0118.6-4631 


289.226 

-69.872 

0.134 

4.5 

1.78 

F 







J0124.6-2322 


188.135 

-81.611 

0.113 

7.5 

2.31 

F 


J0123.7-2312 

0.217 

1.225 


S,F 

J0128.0+6330 


126.998 

0.922 

0.255 

5.0 

2.57 

F 


J0128.4+6257 

0.432 

1.262 


F 

J0129.4+2618 


133.451 

-35.784 

0.333 

4.9 

2.56 

F 

J0127.9+2551 


0.813 

0.707 

T 


J0158.6+8558 


124.201 

23.262 

0.180 

4.2 

2.52 

F 

J0145.6+8600 


0.334 

0.691 

T 


J0214.5+6251C 


132.251 

1.495 

0.134 

4.1 

2.26 

F 






c,F 

J0218.7+6208C 


132.937 

0.975 

0.078 

10.5 

2.77 

F 

J0217.3+6209 


0.208 

0.753 

T 

S,c,F 

J0219.1-1725 

IRXS J021905.8-172503 

191.883 

-67.564 

0.148 

4.3 

1.92 

F 





T 


J0221.3+6025C 


133.810 

-0.528 

0.110 

4.6 

2.45 

F 


J0221.1+6059 

0.391 

1.451 


c,F 

J0221.4+6257C 


132.962 

1.856 

0.118 

9.5 

2.64 

F 


J0223.5+6313 

0.316 

1.11 

T 

^,c,F 

J0227.2+6029C 


134.471 

-0.220 

0.097 

6.1 

2.38 

F 








^All the values reported in these columns are from the 2FGL catalog. 

'^Name of the 3FGL source associated with the 2FGL source with positional coincidence evaluated using dgg.g. 

^Closest 3FGL source having a distance ^ 99.9 < A < 1° from the position of the 2FGL source. 

^In this column is reported the angular separation (A) between the 2FGL source and the 3FGL sources associated using ^ 99.9 or the closest 3FGL source. 

®T: The 2FGL source and one of the initial seeds for the 3FGL analysis have angular separation < dgg,g. 

^S: The 2FGL source was split/resolved into one or more seeds; c: The 2FGL source was flagged with c, i.e., possibly contaminated by the diffuse emission; F: the 2FGL source 
had analysis flags; E: The 2FGL source has a distance < l°from an extended 3FGL source. 

This table is published in its entirety in the electronic edition of the Astrophysical Journal Supplements. A portion is shown here for guidance regarding its form and content. 
Similar tables are available, only in the electronic edition, for lost OFGL, lost IFGL, and lost IFHL sources. 
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list, the catalog construction incorporates most known clouds into the diffuse model, 
and so no sources of this type are identified in the catalog. 


3. Correlated Variability. Variable sources, primarily AGN, whose 7 -ray variations can be 
matched to variability seen at one or more other wavelengths, are considered to be firm 
identifications. Although some cases are well documented, such correlated variability 
is not always easily defined. We conservatively require data in more than two energy 
bands for comparison. Finding a blazar to have a high X-ray flux at the same time as a 
7 -ray flare, for example, does not qualify if there is no long-term history for the X-ray 
emission. We include those sources whose variability properties are documented either 
in papers or with Astronomer’s Telegram^. This list does not represent the result of 
a systematic study. Ongoing work will undoubtedl y enlarge this list. The one Galactic 
source identified in this way is nova V407 Cygni flAbdo et al.ll 2010 hll . Similarly short 
duration tangent gamma-ray emission observed from the classical novae, V959 Mon 
2012 a nd V1324 Sco 2012, were not detected in the 4-year integrated analysis in the 
3FGL l Ackermann et ah 2014a). 


We include one exception to these rules. The Grab PWN is listed as a firm identification 
even though it does not meet any of these criteria. The well-defined energy spectrum, distinct 
from the Grab pulsar spectrum and matc hing spectra seen a t both lower and higher energies 
provides a unique form of identification flAbdo et al.ll 2010 ell . 


In total, we firmly identify 238 out of the 3033 3FGL sources. Among those, 143 are 
pulsars, 66 are active galactic nuclei (BGU, BLL, FSRQ, NLSYl, or RDG) 12 are SNR, 4 
are binaries (BIN or HMB), 9 are PWN, 2 are normal galaxies, 1 is a massive star-forming 
region, and 1 is a nova (Table | 6 ]). 


5.2. Automated Source Associations 


Our approach for automated source association closely f ollows that used for the 2FGL 


and details of the method are provided in lAbdo et al.l fcOlOdl. IFGL) and iNolan et al.l (12012 
2FGL). 


In summary, we use a Bayesian approach that trades the positional coincidence of pos¬ 
sible counterparts with 3FGL sources against the expected number of chance coincidences 
to estimate the probability that a specific counterpart association is indeed real (i.e., a phys- 


^^See http://www.astronomerstelegram.org. 
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ical association). As for IFGL and 2FGL, we retain counterparts as associations if they 
reach a posterior probability of at least 80%. We apply this method to a set of counterpart 
catalogs for which we calibrate the prior source association probabilities using Monte Garlo 
simulations of fake 3FGL catalogs. 

Table [12] lists the catalogs used in the automatic association procedure, organized into 
four categories: (1) known or plausible y-ray-emitting source classes, (2) surveys at other 
frequencies, (3) GeV sources, and (4) identified y-ray sources. The first category allows us 
to assign 3FGL sources to object classes, while the second category reveals multiwavelength 
counterparts that may suggest the possible nature of the associated 3FGL sources. The third 
category allows assessment of former GeV detections of 3FGL sources, and the fourth keeps 
track of all hrm identifications (cf. 15.ip . For this last category we claim associations based 
on the spatial overlap of the true counterpart position with the 3FGL 99.9% conhdence error 
ellipse. 

With respect to 2FGL, we updated all catalogs for which more comprehensive compila¬ 
tions became available. As for IFGL and 2FGL, we separately consider energetic and nearby 
pulsars, with spin-down energy flux E/d? > 5 x 10^^ erg kpc“^ s“^. We also consider mil¬ 
lisecond pulsars (MSPs) separately. For spin period P (s) and spin-down rate P, we define 
MSPs as pulsars satisfying log^g ^ + 19-5 + 2-5 logj^g -P < 0. 

Gatalogs indicated with an asterisk in Table IT^ have source location uncertainties greater 
than those for the 3FGL sources. Gatalogs indicated with a dagger have extended sources, 
with sizes greater than the source location uncertainty regions for 3FGL. For these cata¬ 
logs we cannot apply the Bayesian association method. For the former catalogs we base 
associations on overlap of the 95% confidence error radii. For the latter, we require overlap 
between the extended source and the 95% confidence radius (semi-major axis) of the 3FGL 
source. These approaches are much less reliable than the Bayesian associations, so we do 
not claim any associations based on overlap in our hnal catalog. We record, however, any 
spatial overlap with a TeV source in the FITS file version of the catalog, and use a special 
flag in our catalog (TEVCAT_FLAG), distinguishing point-like (P, angular diameter < 20') from 
extended (E) TeV counterparts (see App. [B|). We furthermore list all unassociated 3FGL 
sources that are spatially overlapping with SNRs or PWNe in Table [T] 
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Table 12. 


Catalogs Used for the Automatic Source Association Methods 


Name 

Objects'^ 

Ref. 

High EjS pulsars 

213 

Manchester et al. 12005')*’ 

Other normal pulsars 

1657 

Manchester et al. (2005)^ 

Millisecond pulsars 

137 

Manchester et al. (2005)^ 

Pulsar wind nebulae 

69 

Collaboration internal 

High-mass X-ray binaries 

114 

Liu et al. (2006^ 

Low-mass X-ray binaries 

187 

Liu et al. (2007^ 

Point-like SNR 

157 

Green (2009) 

Extended SNR^ 

274 

Green (2009) 

0 stars 

378 

Maiz-Aoellaniz et al. (2004') 

WR stars 

226 

van der Hucht (2001) 

LBV stars 

35 

Clark et al. (2005^ 

Open clusters 

2140 

Dias et al. (2002) 

Globular clusters 

160 

^larm (1996) 

Dwarf galaxies^ 

100 

McConnachie (2012) 

Nearby galaxies 

276 

Schmidt et al. fl993^ 

IRAS bright galaxies 

82 

Sanders et al. ('2003) 

BZCAT (Blazars) 

3060 

Massaro et al. ('2009) 

BL Lac 

1371 

V&on-Cettv & Veron (2010) 

AGN 

10066 

V&on-Cettv & V&on (2010) 

QSO 

129,853 

V&on-Cettv & V&on (2010) 

Seyfert galaxies 

27651 

V&on-Cettv & V&on (2010) 

Radio loud Seyfert galaxies 

29 

Collaboration internal 

IWHSP 

1000 

Arsioli et al. f2014) 

WISE blazar catalog 

7855 

D’Abrusco et al. (2014) 

NRAO VLA Sky Survey (NVSS)'' 

1,773,484 

Condon et al. ('1998) 

Sydney University Molonglo Sky Survey (SUMSS)'^ 

211,050 

Mauch et al. ('2003) 

Parkes-MIT-NRAO survey‘d 

23277 

Griffith & Wrieht (1993) 

CGRaBS 

1625 

Jleale^^^al^ ('2008) 

CRATES 

11499 

^leale^^^al^ ('2007) 

VLBA Calibrator Source List 

5776 

http: / / www.vlba.nrao.edu / astro / calib / vlbaCalib.txt 

ATCA 20 GHz southern sky survey 

5890 

^Rirghv et al. ('2010) 

ATCA follow up of 2FGL unassociated sources 

424 

Petrov et al. ('2013) 

ROSAT All Sky Survey (RASS) Bright and Faint Source Catalogs^ 

124,735 

Voges_et al. (1999).'* 

58 months BAT catalog 

1092 

^aumgart^r et al. ('2010) 

^th jgjg catalog 

723 

Bird et al. ('2010) 

1st AGILE catalog* 

47 

Pittori et al. f2009) 

3rd EGRET catalog* 

271 

Hartman et al. fl999) 

EGR catalog* 

189 

Casandiian Sz Grenier ('2008) 

OFGL list* 

205 

Abdo et al. f2009d. OFGL) 

IFGL catalog* 

1451 

Abdo et al. f2010d. IFGL) 

2FGL catalog* 

1873 

Nolan et al. (2012. 2FGL) 

IFHL catalog* 

514 

Ackermann et al. f2013a. IFHL) 

TeV point-like source catalog* 

82 

http: //tevcat.uchicago.edu/ 

TeV extended source catalog^ 

66 

http: //tevcat.uchicago.edu/ 

LAT pulsars 

147 

Collaboration internal 

LAT identified 

137 

Collaboration internal 
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5.2.1. Active Galactic Nuclei associations 

AGN, and in particnlar blazars, are the largest class of associated sonrces in 3FGL at 
high Galactic latitndes. 


In SLAG, as in the Second LAT Gatalog of AGN (2LAG, lAckermann et al.l boilb ). 
we added another association method to the antomatic one described above. This is the 


radio, infrared and optical sources fe.e.. 

de Ruiter et ah 

1977: 

Prestaee & Peacock 

1983: 

Sutherland & Saunders 

1992; 

Lonsdale et al.lll998: 

Masci et ahlboOl: 

Ackermann et ah 

2ni1bi 


It is based on nniform snrveys in the radio and in X-ray bands, enabling us to search for 
possible counterparts among the faint radio and X-ray sources. The LR method makes use 
of counterpart densities (assumed spatially constant over the survey region) through the 
log N — log S relation and therefore the source flux. This approach has been already used 
in 2LAG and we refer the reader to § 3.2 of that paper for a comprehensive description 
of the method, which computes the probability that a suggested association is the ‘true’ 
counterpart @ 

A source is considered as a likely counterpart of the y-ray source if its reliability (see 
Eq. 4 in 2LAG) is greater than 0.8 in at least one survey. 

In total, our automatic association procedure based on the Bayesian method hnds 1663 
3FGL sources that are associated with AGN while the LR-based association method hnds 
1340 3FGL sources. For 405 sources only the Bayesian method provides an association, and 
for 82 sources only the LR method does so. 

Overall, 3FGL includes 1745 sources associated with AGN (58% of all 3FGL sources) 
of which 1145 are blazars, 573 are candidate blazars, 15 are radio galaxies, 5 are Seyfert 
galaxies, and 3 are other AGN. The Seyfert galaxies are n arrow-line Seyfert 1 galaxies that 
have been established as a new class of y-ray active AGN (lAbdo et ahl 1200911) . 


Gomparing to 2FGL we can make the following observations: 


The 3FGL includes 610 more sources of AGN type than the 2FGL, i.e., a 76% increase. 


^®We note that the implementation of the LR method for the 2LAC associations was plagued with an 
error in the management of the sky coordinate precession that affected some of the associations. These 
false associations were also included in the active galaxy associations in 2FGL. The error has been fixed 
and the 2FGL associations re-derived. The corrected 2FGL catalog file has been delivered to the FSSG 
for distribution. In the present work, comparisons with 2FGL findings are based on this corrected set of 
associations. 




































69 


The fraction of new sources (not present in 2FGL) is slightly higher for hard-spectrum 
(i.e., r < 2.2) sources than for soft-spectrum ones (i.e., T > 2.2), 51% vs. 47% 
respectively, but the relative increase reaches 72% for very hard-spectrum (i.e., T < 
1.8) sources. In the 3FGL, 477 counterparts are new (81 FSRQs, 146 BL Lacs, 240 
candidate blazars of unknown type, 10 non-blazar objects); the other counterparts were 
present in previous AGN Fermi catalogs but not included in any of the OFGL, IFGL, 
or 2FGL catalogs for various reasons (e.g., the corresponding y-ray sources were not 
associated with AGN, had more than one counterpart or were flagged in the analysis). 

• The fraction of counterparts of unknown type (named ‘bcu’) has increased notably 
between the two catalogs (from 20% to 28%). The number of these sources in the 
3FGL has increased by more than a factor of 2.5 relative to that in the 2FGL, becoming 
almost equal to that of FSRQs. This increase is mainly due to the lower probability of 
having a published high-quality spectrum available for these fainter sources because of 
the lack of optical observing programs. In 3LAG, sources of the ‘bcu’ type have been 
divided into three sub-types depending of the multi-wavelength information available 
to characterize their ‘blazarness’. In this paper we do not propagate this sub-division 
and refer to the 3LAG for census. 

• The relative increase in ‘bcu’ drives a drop in the proportions of FSRQs and BL Lacs, 
which only represent 29% and 41% of the 3FGL respectively (38% and 48% for 2FGL). 
The relative increase in the number of sources with respect to 2FGL is 34% and 42% 
for FSRQs and BL Lacs respectively. 

• Out of 825 AGN in the 2FGL, a total of 68 are missing in the full 3LAG sample, most 
of them due to variability effects. A few others are present in 3FGL but with shifted 
positions, ruling out the association with their former counterparts. 


5.2.2. Normal galaxies 

The y-ray emission of normal galaxies is powered by cosmic-ray interactions with inter¬ 
stellar gas and radiation. They are numerous but typically faint relative to active galaxies. 
The most luminous of the normal galaxies are starburst galaxies, which have very high den¬ 
sities of gas and massive star formation near their centers. Less distant are normal galaxies 
in the local group. As described above we searched for associations with sources in catalogs 
of nearby galaxies and IRAS bright galaxies (Table [12]). 

In the 3FGL catalog we do not hnd additional associations with normal galaxies rel¬ 
ative to those reported already in 2FGL : starburst galaxies M82 (3FGL J0955.4-1-6940), 
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NGC 253 (3FGL J0047.5-2516) NGG 1068 (3FGL J0242.7-0001), and NGG 4945 (3FGL 
J1305.4—4926), and local gronp galaxies LMG (3FGL J0526.6—6825e), M31 (3FGL J0042.5+4117), 
and the SMG (3FGL J0059.0-7242e). 

Five sonrces in the 3FGL catalog lie within the extended sonrce model for the LMG and 
are otherwise nnassociated with connterparts at other wavelengths. These sonrces (3FGL 
J0456.2-6924, 3FGL J0524.5-6937, 3FGL J0525.2-6614, 3FGL 0535.3-6559, and 3FGL 
J0537.0—7113) are classified here as ‘gah based solely on the spatial coincidence with the 
LMG and their associations are listed as LMC field. Their particnlar natnres remain nncer- 
tain. 


5.2.3. Pulsars 


Becanse pnlsed emission can be snch a clear signatnre, pnlsars represent the largest 
class of firmly identified astrophysical objects in the 3FGL catalog. An extensive discnssion 
of 7 -ray pnlsar prop erties is fonnd in the Second Fermi Large Area Telescope Gatalog of 
Gamma-ray Pnlsars jAbdo et al.ll2013l. 2PG). The pnblic catalog of LAT-detected pnlsars 
is regnlarly npdatecJ^ At the time of the 3FGL association analysis, this catalog had 
147 pnlsars (Table [12]). Only 137 of the LAT-detected pnlsars have associations in the 
3FGL catalog, however (Table Ej). The missing ten did not reach the TS > 25 criterio n 
based on their average flnxes. Three of these are PSR J0737—3039A flGnillemot et al.ll2013l) . 
J1640+2224, and J1705—1906 (IHon et al.ll20141) . and the remaining seven are flagged in the 
2PG ‘spectral resnlts’ tables as being either too faint to fit, or reqniring an on-peak analysis 
to obtain spectra. 


5.2.4- Pulsar wind nebulae 

In addition to the fonr pnlsar wind nebnlae fonnd in 2FGL (Grab, Vela-X, MSH 15—52, 
HESS J1825—137), the 3FGL catalog inclndes seven new PWNe associations. Five of these 
are firm identifications becanse they are spatially extended LAT sonrces (see Table dj): 
HESS J1303-631 (3FGL J1303.0-6312e), HESS J1616-508 (3FGL J1616.2-5054e), HESS 
J1632-478 (3FGL J1633.0-4746e), HESS J1837-069 (3FGL J1836.5-0655e), and HESS 
J1841—055 (3FGL J1840.9—0532e). The other two are positional associations with known 
PWNe: G279.8-35.8 (3FGL J0454.6-6825) and GO.13-0.11 (3FGL J1746.3-2859c). 


^®See https;//confluence.slac.Stanford.edu/display/GLAMCDG/Public+List+of+LAT-Detected+Gamma-Ray+Pulsa] 
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5.2.5. Globular clusters 

Two globular cluster associations from the 2FGL catalog are not found in the 3FGL 
catalog: 


1. 2FGL J1727.1—0704, previously associated with IG 1257, is found as 3FGL J1727.6—0654. 
This source is not formally associated with the globular cluster. 

2. 2FGL J1824.8—2449, which was associated with NGG 6626, has been hrmly identified 
as PSR J1824—2452A f .Tohusou et ah 2013 ). Its catalog listing is 3FGL J1824.6—2451. 


The number of globular clusters associated with LAT sources does continue to grow. 
New associations are NGG 2608 (3FGL J0912.2-6452), NGG 6316 (3FGL J1716.6-2812), 
NGG 6441 (3FGL J1750.2-3704), NGG 6541 (3FGL J1807.5-4343), NGG 6717 (3FGL 
J1855.1—2243), and NGG 6752 (3FGL J1910.7— 6000). NGG 6752 had previously been 
noted as a likely LAT source by Tam et ah f 2011bl) . 


5.2.6. Supernova remnants 


Twelve SNRs are firmly identified in the 3FGL catalog as spatially extended sources 
(see Tabled]). Six had previously appeared in the 2FGL catalog: IG 443, W 28, W30, W44, 


W51 G, and the Gygnus Loop. Additio ns are: S147 (3FGL J0540.3+2756e, iKatsuta et ah 


2 OI 2 II. PuDDis A Id FGL J0822.6-4250e, iHewitt et al.ll2012l L Vel a Jr. f3FGL J0852 .7-4631e. 

Tanaka et ah 2011 ). RX J171 3.7—3946 (3FGL J 1713.5—3945e, Abdo et ah 2011d ). Gamma _ 

Gygni (3FGL J2021.0+4 031e, Lande et aLll2012 '). and HB21 (3FGL J2045.2+5026e, Reichardt et af 
2012 : Pivato et ah 2013). 


Additionally we consider eleven unresolved 3FGL sources as being confidently associated 
with SNRs, based on individual studies of these SNRs in LAT data fsee lFerrand fc Safi-Harb 
201^. and references therein)@. These are given the ‘snr’ designator (Table Ej). The 2FGL 
sources that were designated ‘snr’ have been added to the ‘spp’ class in 3FGL. Many of 
the 68 SNR or PWNe in this table are spatially extended sources at other wavelengths, 
and therefore the chance probability of an overlap with a LAT source is non-negligible. As 
for previous LAT catalogs, we encourage great care in any analysis using these potential 
associations. 


http: //WWW. physics . umanitoba. ca/snr/SNRcat/ 
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5.2.1. Binaries 


Three HMB sources that appeared in 2FGL are also found in 3FGL; LS I +61 303 
(3FGL J0240.5+6113), IFGL J1018.6-5856 (3FGL J1018.9-5856), and LS 5039 (3FGL 
J1826.2—1450). All were hrmly identihed by binary periodicity. We note that each of the 
three is a TeV emitter but that the two other binary systems that have been detected in 
the TeV energy regime, HESS J0632+057 and PSR B1259—63, do not have counterparts in 
3FGL. PSR B1259—63/LS 2883 is a bright LAT source during a small part of the 3.4-year 
binary period following periastron ( Tam et ah 2011a : Abdo et ah 2011a). HESS J0632+057 
has not been detected at all. 

A fourth HMB from 2FGL, Gygnus X—3 (2FGL J2032.1+4049) does not appear in 
3FGL. It is an intermittent LAT source and was not active enough averaged over the four 
years of this catalog to produce a signihcant detection. Eta Garinae, which appeared in 
2FGL as a possible massive star association, is included as 3FGL J1045.1—5 941, identihed as 


a binary sy stem. The LAT data exhibit the known 5.5-year binary period flReitberger et ah 


2 OI 2 I. I 2 OI 41 I . Although liste d in the catalog as a nova, V407 Gygni is also a binary system. 
Its flaring activity in 2010 (lAbdo et aLll2010hll was bright enough that it appears as 3FGL 
J2102.3+4547. 


5.2.8. Multiwavelength associations 


In addition to the catalogs of classified sources, we also searched for associations with 
catalogs of radio and TeV sources. Our association procedure for AGN heavily relies on 
associations with radio sources as most of the 7 -ray emitting AGN are bright sources of 
radio emission (see § I5.2.ip . In fact, essentially all of the radio associations we find have 
been classified subsequently as AGN. 


We did not search for general associations with infrared, optical, or soft X-ray catalogs. 
Within the LAT source error regions we would find multiple potential counterparts, most of 
which necessarily would be due to chance, since many of the sources in these catalogs are ther¬ 
mal in nature. We included, however, the hard X-ray catalogs INTEGRA L-IBIS and Swift- 


BAT, and blazar candidates extracted from the infrared WISE catalog fjP’Abrusco et ah 


20141 1 in the automated association pipeline. These data, when included in a study of the 


spectral energy distributions to evaluate their synchrotron peak frequencies and general be¬ 
haviors, help in understanding the natures of the candidate counterparts. This was done 
especially for all the sources classified as bcu/BGU and agn. 
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5.2.9. Statistics of association results 

In total we find that 1976 of the 3033 sources in the 3FGL catalog (59%) have been 
associated with at least one non-GeV 7 -ray counterpart by the automated procedures. Table 
[13] summarizes the association results. 


5.3. GeV and TeV Source Associations 


Through 2014 August, 155 flaring Fermi-LAT sources were detected and promptly re¬ 
ported in more than 249 Astronomer’s Telegram^. Of these sources, six are not in 3FGL. 
For two of these the flaring state was detected outside the time interval covered by 3FGL: 
S5 1044-1-71 (a 2FGL source classihed as an FSRQ) and PMN J1626—2426 (in the prox¬ 
imity of an unassociated 2FGL source flagged as potentially contaminated by the diffuse 
emission). The other four are Gyg X—3 (an HMB associated with 2FGL J2032.1-1-4049; 
see § I5.2.7p . GGRaBS J1848-I-3219 (an FSRQ associated with 2FGL J1848.6-1-3241), PKS 
1124—186 (an FSRQ associated with 2FGL J1126.6—1856) and PKS 2123—463 (an FSRQ 
associated with 2FGL J2125.0—4632). The reason that these three FSRQs are missing from 
the 3FGL catalog is probably that they have average fluxes below the detection threshold. 


Sources in 3FGL that are positionally associated with sources seen by the ground-based 
TeV telescopes are of particular interest for broad-band spectral studies. As for the 2FGL 
catalog we studied associations with the TeVGat@ compilation of sources detected by very- 
high-energy observatories. The energy threshold for TeVGat sources is not uniform, but it 
is typically greater than 100 GeV. We used a compilation of TeVGat sources prepared on 
2014 October 27 that has 148 unique entries. This comprises the so-called “Default” and 
“Newly Announced” TeV catalogs. We note that, as in IFGL and 2FGL, the “Galactic 
Gentre Ridge” was not included for association pu rposes since it represen ts diffuse emission 
over an extended region along the Galactic plane flAharonian et al.ll2006h . 


This TeVGat compilation by its nature does not represent a complete survey, and our 
general statistical procedure for evaluating probabilities of chance association could not be 
applied. As for 2FGL we separately considered point-like and extended TeVGat sources 
(Table [T2|). For point-like sources the criterion for association was overlapping 95% source 
location regions (indicated by ‘P’ in the ‘TeV’ column of Table Hj). For extended sources 
(LAT or TeVGat) the criterion was spatial overlap within their respective angular extents 


^^https://www-glast.Stanford.edu/cgi-bin/pub_rapid 
^^http://tevcat.uchicago.edu 
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(indicated by ‘E’). We note that, in the literature, the shapes of the extended TeV sources 
are usually approximated to a circle or to an ellipse. For the purposes of our association 
pipeline, we imposed a circular geometry on all extended TeV sources, setting the radius to 
the length of the semi-major axis. In the case of TeV sources whose morphologies depart 
significantly from a simple ellipse or a circle, this simplification of their geometry could lead 
to missed 3FGL associations. 


In total, 124 3FGL sources have TeV counterparts. Of the 148 TeV sources considered, 
117 have 3FGL associations and, out of these, 6 TeV sources have multiple 3FGL associations. 
Five of these are extended TeV sources and one is the Grab, which is associated with both 
the synchrotron and the inverse-Gompton Grab 3FGL sources (3FGL J0534.5-|-2201s/3FGL 
J0534.5-|-2201i). The TeV sources HESS J1018—589, Westerlund 2, HESS J1458—608 and 
MGRO J2019-I-37 have two 3FGL associations each while Westerlund 1 has three 3FGL 
associations. 


We note that the TeV source HESS J1018—589 has two components, denoted A (a point 
source) and B (extended emission). Of the two 3FGL associations, 3FGL J1018.9—5856, a 
LAT high-mass binary, lies closer to location A and 3FGL J1016.3—5858, a LAT pulsar, lies 
closer to location B. 


Table [H] shows the associations between extended TeVGat sources and 3FGL catalog 
sources. Some of these, designated with e appended to their source names, were explicitly 
modeled as extended sources corresponding to H.E.S.S. sources (see li 13.41) . 

Out of the 58 TeV AGN, 57 have associated sources in 3FG L. Only HESS J1943+213 


tentatively cla.ssified as a high-synchrotron peaked (HSP) blazar flHESS Gollaborationll2011 


Leahy fc TianI 120121) . does not have a 3FGL association. This presumed blazar is unique 
in the TeV sky in that it shines through the Galactic plane. We note that the VERITAS 
source VER J2016-I-371 is positionally associated with 3FGL J2015.6-1-3709 although the 
VERITAS source is probably a PWN flAliu et al.l 120141) and the LAT source is associated 
with an FSRQ of unknown redshift. 

The Milagro source, MGRO J2031-I-41, which comprises TeV J2032-I-4130, was also in¬ 
cluded as a separate source in the TeV list that was used to evaluate the assoc iations with 


3FGL because it is postulated that its emission is due to m ore than one source flAbdo et ah 


20071) . Due to the large extent of MGRO J2031-I-41 (1?8, lAbdo et al.ll2012al) . it has posi¬ 
tional overlap with eight 3FGL sources in addition to 3FGL J2032.2-1-4126. We have listed 
just one LAT source 3FGL J2028.6-|-4110e, the Gygnus Gocoon, as being associated with 
MGRO J2031+41. 


Due to its large extent (2?6), the Milagro Geminga source flAbdo et al.l l2009hl) has 
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positional overlap with two 3FGL sources, 3FGL J0633.9+1746 (the Geminga pulsar) and 
the unidentihed source 3 FGL J0626.8+1743. We have associated only the Geminga pulsar 
with the Milagro source. lAbdo et ahl fj2009hl) postulate that the Milagro emission could be 
due to a pulsar-driven wind associated with Geminga. 


The TeV sources Boomerang and SNR G103.3+02.7 have positional overlap at TeV en¬ 
ergies. They are each positionally coincident with the same two 3FGL sources, the uniden¬ 
tified source 3FGL J2225.8+6045 and the pulsar, 3FGL J2229.0+6114. We have associated 
SNR G103.3+02.7 with 3FGL J2225.8+6045 and Boomerang, classified as a PWN at TeV 
energies, with 3FGL J2229.0+6114, the LAT pulsar. 

Relative to 2FGL eight TeV sources are newly associated with LAT sources. All but 
one of these sources (HESS J1641—463) had already been detected at TeV energies when 
2FGL and IFHL were released. None of these eight sources, however, have counterparts in 
those catalogs: 


1. The H.E.S.S. Galactic Genter source is associated with 3FGL J1745.6—2859c. The 
corresponding source 2FGL J1745.6—2858, had a large enough position offset that it 
was not considered associated with the TeV source. The Galactic center remains, 
however, a particularly complex region whose detailed study is beyond the scope of 
this paper. 

2. Three HSP blazars (SHBL J001355.9-185406, lES 0229+200 and lES 0347-121) are 
associated with 3FGL sources. 

3. The shell SNR Tycho is associated with the faint LAT source 3FGL J0025.7+6404. 
Although the center of the TeV emission is offset by 0?12 from the LAT source, the 
relatively large uncertainty of the LAT position indicates sufficient positional overlap 
for association. 


The TeV PWN HESS J1809—193 is another new TeV association in 3FGL. Discovered 
at TeV energies in 2007, it is a relatively bright Te V source with an integ ral flux 14% 
that of the Grab Nebula in the same energy band f|Aharonian et ahl 120071) . 


5. The unidentified TeV sources HESS J1626—490 and HESS J1641—463 are also new 
TeV associations for 3FGL. 


Thirty-one TeV sources have no counterparts in 3FGL; one of these, the unidentified 
extended TeV source HESS J1857+026, has an association in IFHL (IFHL J1856.9+0252). 
It is the only TeV source to have a IFHL association but none in 2FGL or 3FGL. The other 
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TeV sources having no counterparts in 3FGL comprise 10 other unidentihed TeV sources, 
10 PWNe, 2 binaries, 4 shell-type SNRs, 2 SNR-molecular cloud associations, 1 composite 
SNR (i.e., 9 out of the 19 SNRs in TeVCat) and the HSP blazars discussed earlier. 

We note that two TeV sources, HESS J1634—472 and SNR W49B, were each associated 
with LAT sources in 2FGL and in IFHL but do not have associations in 3FGL. A source 
coincident with W49B has been detected by the LAT, 3FGL J1910.9-1-0906 but this does 
not have positional overlap with the TeV detection of W49B. 

The TeV source, HESS J1427—608 had an association in 2FGL but was not associated 
with a LAT source in IFHL and does not have an association in 3FGL. 



Table 12—Continued 


Name Objects^ Ref. 


^Number of objects in 
the catalog. 

'^http: / / WWW. atnf. csiro.au / research / pulsar / psrcat 

"^All-sky surveys used 
only in the Likelihood 
Ratio method, see ^ 15.2.11 

^The RASS Faint 
Source Catalog 

is available from 

http://www.xray.mpe.mpg.de/rosat/survey/rass-fsc/ 








Table 13. Statistics of Source Associations 


Gategory 

OFGL 

IFGL 

2FGL 

IFHL'* 

3FGL 

Total 

205 

1451 

1873 

514 

3033 

Associated 

168 

821 

1224 

449 

2023 

Unassociated 

37 

630 

649 

65 

1010 

New 7 -ray sources^ 

- 

1265 

762 

52 

1312 

Sources associated with former LAT detections 

- 

186 

1111 

462 

1721 

Sources associated with former GeV detections^ 

74 

162 

170 

4 

206 

Firmly identified sources 

31 

65 

124 

60 

238 

Sources associated with at least one object of known type 

153 

623 

952 

385 

1398 

Sources which have counterparts only in the multiwavelength catalogs 

- 

92 

214 

58 

576 


=-lFHL: >10 GeV 

^Non-overlapping 95% source location confidence contours compared to previous LAT catalogs, at the level of overlapping 
95% source location confidence contours. 

^Here only the lAGL, 3EG, and EGR catalogs are considered. 
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Table 14. Associations of 3FGL with Extended TeV Sources 


TeVCat Name'^ 3FGL Name 


Boomerang 
CTA 1 
CTB 37A 
Geminga 
HESS J1018-589 
HESS J1303-631 
HESS J1356-645 
HESS J1458-608 
HESS J1503-582 
HESS J1507-622 
HESS J1614-518 
HESS J1616-508 
HESS J1626-490 
HESS J1632-478 
HESS J1640-465 
HESS J1708-443 
HESS J1718-385 
HESS J1745-303 
HESS J1800-240 
HESS J1804-216 
HESS J1809-193 
HESS J1825-137 
HESS J1834-087 
HESS J1837-069 
HESS J1841-055 
HESS J1848-018 
HESS J1858+020 
IG 443 

Kookaburra (Rabbit) 
Kookaburra PWN 
MGRO J1908+06 
MGRO J2019+37 
MGRO J2031+41 
MSH 15-52 
RX J0852.0-4622 
RX J1713.7-3946 
SNR G292.2-00.5 
Terzan 5 
TeV J2032+4130 
Vela X 

VER J2019+407 
W 28 
W 51 

Westerlund 1 


Westerlund 2 


J2225.8+6045, J2229.0+6114 
J0007.0+7302 
J1714.5-3832 
J0633.9+1746 

J1016.3-5858, J1018.9-5856 

J1303.0-6312e 

J1356.6-6428 

J1456.7-6046, J1459.4-6053 

J1503.5-5801 

J1506.6-6219 

J1615.3-5146e 

J1616.2-5054e 

J1626.2-4911 

J1633.0-4746e 

J1640.4-4634C 

J1709.7-4429 

J1718.1-3825 

J1745.1-3011 

J1800.8-2402 

J1805.6-2136e 

J1810.1-1910 

J1824.5-1351e 

J1834.5-0841 

J1836.5-0655e 

J1840.9-0532e 

J1848.4-0141 

J1857.9+0210 

J0617.2+2234e 

J1418.6-6058 

J1420.0-6048 

J1907.9+0602 

J2021.1+3651, J2017.9+3627 

J2028.6+4110e 

J1514.0-5915e 

J0852.7-4631e 

J1713.5-3945e 

J1119.1-6127 

J1748.0-2447 

J2032.2+4126 

J0833.1-4511e 

J2021.0+4031e 

J1801.3-2326e 

J1923.2+1408e 

J1648.3-4611, J1650.3-4600, 
J1651.5-4626 

J1023.1-5745, J1024.3-5757 
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5.4. Properties of Unassociated Sources 


Among the 3033 sources in the 3FGL catalog, 2033 have associations or identifications 
with known astrophysical objects. Although that number is greater than the total number of 
sources in the 2FGL catalog, 1010 (33%) of the 3FGL sources remain unassociated. Among 
these unassociated sources are many that were found in previous LAT catalogs, indicating 
that some persistent mysteries remain despite extensive efforts to find associations over the 
past few years. The continued prevalence of unassociated sources is expected, as the improve¬ 
ment in sensitivity with four years of flight data and improvements to the characterization 
of backgrounds have allowed Fermz-LAT to probe the y-ray sky to unprecedented depths. 
As a result, direct comparison to previous releases is difficult. 


The distribution of unassociated sources on the sky is compared in Figure |23] to the dis¬ 
tribution of the associated sources. The plot reveals some important features that should be 
kept in mind when considering unassociated 3FGL sources. Of the 992 unassociated sources 
in the 3FGL, 334 fall within the Galactic plane (|6| < 5°). This leaves 658 unassociated 
at |6| > 5° sources for an average density of 1.75 x 10“^ sources deg“^. While this density 
reflects the finite angular resolution and sensitivity of the instrument, some of the shortfall 
of associations is attributable to the fact that not all areas of the sky have been mapped 
uniformly at other wavelengths. For instance, a 3269 deg^ “overlap region” in the North 
Galactic Gap intensively covered by five radio catalogs GB6 (6 cm) FIRST (20 cm), NVSS 
(20 cm), WNSS (92 cr n), and VLSSr (4 m), as well as by the Sloan Digital Sky Survey (SDSS 
DR9) optical survey flKimball &: Ivezid 120081) . contains only 31 unassociated sources. The 
corresponding density of unassociated sources is 9.48 x 10“^ sources deg“^, nearly half of the 
overall average. 


Within the Galactic plane, the unassociated source population is a combination of both 
Galactic and extragalactic source classes. Outside the plane, the LAT-detected AGN source 
density is 0.045 deg“^. By extrapolation, this implies that there should be 160 detectable 
AGN within the 10-degree band surrounding the Galactic plane (not accounting for incom¬ 
pleteness of AGN catalogs at low latitudes). Only 56 sources are associated with active 
galaxies in this region. At low latitudes the LAT detection threshold is higher, and catalogs 


Table 14—Gontinued 


TeVCat Name^ 3FGL Name 


^From 

http://tevcat.uchicago.edu 
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Fig. 23.— Distributions in Galactic latitude b of unassociated sources (shaded red region), all 
associated sources (blue histogram) and all active galaxy source classes (black line). Binned 
in sin 6, an isotropic distribution would be flat. 

of active galaxies are incomplete, but extrapolation from higher latitudes suggests that fewer 
than ~100 of the unassociated sources in the region |6| < 5° could be active galaxies. 

Of the remaining unassociated sources, we expect most to be Galactic objects. If we 
use the fractions of Galactic associated sources in this latitude range as a guide, we hnd that 
nearly half should be pulsars (47%), and most of the remainder should be SNRs (44%), with 
only a small number of other Galactic sources (9%). These yet-to-be-detected pulsars may 
never be seen to pulse: part of the neutron star population will have magnetic and rota- 
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tion a xis orientations c ausing 7 -ray emission to be spread over a large part of the rotation 
phase flHou et alJl2014j) . This, combined with a low 7 -ray flux in the presence of high back - 
ground, renders blind periodicity searches insensitive to pulsations flDormodv et alJl201lh . 
That many of these sour ces are likely to be previously undetected SNRs is also reasonable. 


Ackermann et al.l (l2012bl ) recognized that the distribution of unassociated LAT sources near 


the Galactic plane matched the scale height appropriate for Population I objects, such as 
the SNR parent population of massive stars. 


It should be emphasized that a substantial fraction of the unassociated sources (40%) 
has at least one analysis flag set (§ 13.Op . We hnd that 57% of the sources with at least one 
flag have | 6 | < 5°, which reflects the complexity of the Galactic diffuse emission flj 13.7.31 and 
§ 13.8p . Because of the difficulties in source detection against the bright diffuse background 
along the Galactic plane, unassociated sources with analysis flags set should be considered 
with caution. That difficulty is acute in the Galactic bulge (within 5° of the Galactic center), 
in the Vela (/ ~ 268°) and Gygnus (/ ~ 80°) regions. The immediate vicinity of the Galactic 
center is particularly uncertain, with 7 sources within 1°, on top of bright diffuse emission. 
We did not attempt in the framework of 3FGL to devise models dedicated to those regions, 
so the source positions and characteristics there are not as reliable as in the extragalactic 
sky. 


We also note a number of clusters of unassociated sources, mostly near the Galactic 
plane. Obvious ones are near (/, 6 ) = (133.5,-1-1), (340,-2). Since there is no specihc search 
for extended sources in 3FGL, many of those clusters are probably extended sources (super¬ 
nova remnants, pulsar wind nebulae, star formation regions). The extended sources which 
are currently declared (§ 13.4p are all relatively bright, with a median signihcance of 32cr. 
This leaves a lot of room for fainter extended sources. 


6. Galactic Source Number Counts 


The 3LAG companion paper discusses briefly the source number counts of extragalactic 
sources. Here we address the Galactic source number counts, following the analysis for 
the IFHL catalog which was based on the method described bv IStrond (120071 1. and which 
addressed energies above 10 GeV. For 3FGL we use energies above 1 GeV. Photon fluxes over 
that energy range (F 35 of Table 0]) are more accurate than over the full band, as explained 
in the 2FGL paper. The much larger number of sources in 3FGL compared to IFHL means 
that the entire analysis is more robust. 


The motivation for performing a Galactic source population analysis is hrstly to obtain 
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estimates of the global source characteristics, i.e., space density and luminosity function, sec¬ 
ondly to estimate the contribution from sources below the detection threshold to the Galactic 
‘diffuse’ emission, and thirdly to generate templates of this emission, to be incorporated into 
diffuse emission models for future source catalogs. Source population analysis also puts the 
detected sources in the context of the total source content of the Galaxy. 

The method is guided by properties of known sources such as pulsars but does not at¬ 
tempt physical modeling of the sources; the approach is essentially geometrical, nevertheless 
the analysis reveals the basic properties of the source population. We refer to the modeling 
of the source population(s) as population synthesis. Since the population synthesis includes 
all sources down to arbitrarily low flux levels (for a given model), it can also be used to study 
the flux limit of the actual catalog, and assess how the observed source number counts are 
affected by the detection procedure. This serves as a consistency check on other methods of 
assessing the detection threshold. 

An essential principle is to use the fact that low-latitude sources probe the high- 
luminosity, low space densities at large distances, while the high-latitude sources constrain 
the low-luminosity, high space density nearby objects. This is because high-luminosity 
sources are rare but visible to large distances which are only sampled in the Galactic plane, 
while low-luminosity sources are common but only visible when nearby, so dominate outside 
the plane. These samples are complementary and allow the full luminosity function to be 
estimated. 


6.1. Source Population Synthesis 


The population synthesis and subsequent analysis is performed using the GALPLOT 

-1 


.y be the luminosity of a source in ph s " m 


software, which is publicly availablJ^. Let L 
some energy range. We use photon luminosities since they are most directly related to the 
detectability of sources and detection thresholds. The luminosity function at Galactocen- 
tric distance R and distance from Galactic plane 2 ; is the space density of sources per unit 
luminosity p{Ly, R, z). The shape of the luminosity function is assu r ned in dependent of po¬ 
sition, i.e., p{Ly,R,z) is separable in Ly and {R,z). After [Strong (1200711 we assume that 


the luminosity function depends on luminosity as for Ly^^in < Ly < Ly^^nax and is zero 
outside these limits. The total space density of sources is p{R, ^) = f p{Ly, R, z) dLy, which 
we normalize to the value p© at {R, z) = (i?©, 0). For a source of luminosity Ly at distance 
d the flux is Sy = Ly/And^. The differential source number counts are dehned as N{Sy) 


: //sourcef orge . net/projects/galplot 
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sources per unit flux over the area of sky considered. At lower S.y, both the luminosity 
function and the spatial boundaries influence N{S^). For this analysis the source fluxes are 
binned in log(S'.^) so that the plotted distributions are proportional to S^N{Sj); we use 5 
bins per decade of S^, appropriate to the statistics available. We use standard Monte Carlo 
techniques to sample p{L^,R,z) throughout the Galaxy, using oversampling to reduce sta¬ 
tistical fluctuations if necessary. We use the sources generated from such simulations to form 
simulated catalogs extending below the 3FGL flux limit and compare the flux distributions 
with the observations. 

In the present work we do not explicitly account for the source detection efficiency in 
3FGL, but simply compare the predictions with the data mindful of the range of the esti¬ 
mated detection threshold. Since the detection efficiency and threshold depend on direction, 
mainly because of the Galactic diffuse emission, and also exposure variations, accounting for 
this would require considerably more study than possible in the present work, and is not 
required for the scope here. 

A large number of models was generated and compared with the data; for this paper we 
choose one reference model which is found to reproduce the data satisfactorily, but which is 
not unique. This suffices to illustrate plausible properties of the Galactic source population. 
A complete study of source number counts, optimizing the model over all parameters and 
considering spatial distributions and more sky regions, or using information about particular 
source classes, is beyond the present paper, and is foreseen in a future work. 


6.2. Model and Comparison With Data 


Our reference model for the luminosity function has p© = 100 kpc“^, and an ® 
dependence on lur ninosity i n the range 2 x 10^^ —2 x 10^® ph s“^ above 1 GeV. The luminosity 
law is discussed in Stronel (2007); the chosen power-law index 1.8 is larger than expected for 
normal pulsars (1.5 or less) or MSPs, but here we wish to encompass both of these source 
types, and also other sources like SNR, with a single power-law function, for simplicity. This 
index is required to fit N{S^) at both low and high latitudes; the steep slope ensures enough 
low-luminosity sources to match high-latitude number counts. 


The distribution in Galactocentric distance is based on the model of Lorimer et ah 


(I 2 OO 6 I) for the distribution of pulsars, taken as representative of Galactic sources. We adopt 
an exponential scale height of 500 pc, guided by that of pulsars; the source number count 
distribution N{S^) depends only weakly on the scale height. This distribution peaks near 
R = 4 kpc and falls to zero at i? = 0; it was chosen for illustration and has not been optimized 
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for the 3FGL source number counts. The spectrum of sources was taken as a power-law with 
exponential cutoff, index 1.4 and cutoff energy 3.2 GeV, with correlated dispersion in these 
parameters as found for pulsars. This is representative of pulsars but is not critical for this 
work since we use integral photon fluxes above 1 GeV as the basis for the analysis of N{Sj). 
The spectrum is used here only when comparing with the spectrum of interstellar emission, 
to estimate the contribution of sources to the diffuse emission (section 6.3). 

Figure 1241 compares the simulated N{S^) with the observed flux distributions of 3FGL 
sources at low latitudes in the inner Galaxy (|6| < 10°, 300° < I < 60°) and high latitudes 
(|6| > 10°), as well as for the full sky for reference. The predictions agree reasonably with 
Galactic plus unassociated sources at low latitudes, and with Galactic associated sources at 
high latitudes (where unassociated sources are probably mainly AGN). At high latitudes it is 
important that the model does not over-predict Galactic plus unassociated sources, and this 
condition is satisfied. The reference model is consistent with the low-latitude source number 
counts, having the observed dependence on flux above the source detection threshold; the 
slope reflects the spatial distribution (independent of the shape of the luminosity function) 
above 10“® ph cm“^ s“^, while the distributions for both the model and observed source 
number counts flatten at lower fluxes, following the luminosity function. 


Some of the unassociated sources in the inner Galaxy may be AGN; using high-latitude 
AGN and unassociated sources (which are probably mainly AGN) and scaling by solid angle, 
we estimate about 30 AGN above the threshold for the inner Galaxy. This can be compared 
to 254 Galactic identihed plus associated sources in the inner Galaxy for the same threshold, 
so that the AGN contribution is negligible, especially considering that most of these should 
be identified/associated AGN and hence excluded from our sample. For comparison there 
are 33 (identified/associated) AGN in this region, similar to our estimate from high latitudes, 
so the selection of Galactic classes certainly avoids the presence of more than a few AGN 
among the unassociated sources. 


_ ^Pulsars, including MSPs, have a range of luminosities 10^^ — 10^” erg s“^ flAbdo et ah 

2013|), corresponding to about 3 x 10^^ — 3 x 10^® ph s“^ (taking a mean energy 2 GeV = 
3.2 X 10“^ erg). This is consistent with the range we have found from our N{S^) analysis at 
low and high latitudes, although we are not assuming anything about the physical nature of 
the sources. 


In this model there are 2.9 x sources in the Galaxy (in the luminosity range con¬ 
sidered), with a total luminosity > 1 GeV of 2 x 10^° ph s“^, or about 6 x 10^^ erg s“^. 
3FGL contains about 266 Galactic sources (identified plus associated in Table E]) so that the 
LAT detects about 1% of the sources in the Galaxy; allowing for a signihcant number of 
unassociated sources at low latitudes being Galactic, a larger number is certainly included 
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(a) Inner Galaxy (b) High Latitudes 




(c) All-Sky 



Fig. 24.— Dependence of source number counts (number of sources per 0.2 dex) on source 
photon flux S above 1 GeV. The markers are source number counts from the 3FGL catalog; 
blue triangles are identihed and associated Galactic sources, red circles are identified and 
associated Galactic, and unassociated sources, and black squares are all sources including 
extragalactic (for reference). The curves are from the reference model described in the text, 
(a) inner Galaxy (|6| < 10°, 300° < I < 60°); (b) high latitudes (|6| > 10°, all longitudes); 
(c) all-sky. 
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in 3FGL. Figure |2l] shows that the distribution of simulated sources (in the reference model) 
continues down to fluxes ~100 times below the detection threshold, the cutoff being due 
to the hnite spatial extent of the Galaxy. The ratio of total flux below threshold to above 
threshold is about 0.25, which gives an estimate of the contribution of the undetected sources 
to the ‘diffuse’ emission (see below). 


Using the work of lWatters fc Romanil fj201ll) we can estimate the number of 7-ray pulsars 
in the Galaxy; they give a pulsar birthrate of 1 per 59 years, which corresponds to 1.7 x 10^ 
pulsars up to age 1 Myr. This is consistent with our model; we include other classes of 
sources, but this shows that our value is plausible. 


6.3. Contribution of Undetected Sources to Diffuse Galactic Emission 


Judging from the turnover in the inner Galaxy N{S^) data, the detection threshold 
there is about 1 x 10“® ph cm“^ s“^, and we adopt this for the following estimates. For the 
inner Galaxy, the total flux from sources is 3 x 10“® ph cm“^ s“^, with 2.4 x 10“® ph cm“^ s“^ 
above threshold, 0.6 x 10“® ph cm“^ s“^ below threshold. So 20% of the total source flux 
is below threshold, and the ratio of flux below/above threshold is 25%. We can use the 
measured diffuse spectrum directly, comparing with the contribution from sources relative 
to interstellar emission: for the inner Galaxy as dehned here, this gives 12% from sources 
above threshold, 3% from sources below threshold, at 1 GeV. These estimates are clearly 
model-dependent, in particular the adopted luminosity function gives a large number of 
low-luminosity sources, but they are certainly of the correct order, since varying the models 
within the range consistent with the data does not change the estimates greatly; details are 
beyond the scope of this paper. 


For comparison with our estimates, IWatters fc Romanil fl201l[) used physical modeling of 
young pulsars, and estimated their contribution to diffuse emission as 2.8%, however using all¬ 
sky averages and for only 6 months of LAT data taking. A study of t he MSP contribution to 
the G alactic emission, for energies above 100 MeV, has been given bv iGregoire fc Knodlseder 
(I2OI3I) : they hnd the contribution is at the few percent level. 


Population synthesis can be used to estimate the increase in the number of sources with 
improved detection limits; in this model, reducing the threshold by a factor 2 would yield 
about twice as many sources at low latitudes. 


Finally we consider the global picture. The luminosity of the interstellar emission from 
cosmic-ray interact ions is about 1 x 10^^ ph s“^ or 4 x 10^® erg s“^ for energies above 1 GeV 
flStrong et al.ll2010j) so that sources have about 20% of the interstellar luminosity. This is also 













an estimate of the contribution of all Galactic sources to the total Galactic 7 -ray intensity, 
averaged over the sky. 


7. Conclusions 

The third Fermi LAT catalog is the deepest-yet in the 100 MeV-300 GeV energy range. 
The increased sensitivity relative to the 2FGL catalog is due to both the longer time interval 
(4 years vs. 2 years for 2FGL) and the use of reprocessed Pass 7 data, which provides a 
narrower PSF above 3 GeV. The 3FGL catalog also benehts from higher-level improvements 
in the analysis, including an improved model for Galactic diffuse emission, and a rehned 
method for source detection. 

The 3FGL catalog includes 3033 sources. The sources are detected (TS > 25) based on 
their average fluxes in the 4-year data set; 647 of the sources are found to be signihcantly 
variable on monthly timescales. We flag 78 (2.6%) of the sources as potentially being related 
to imperfections in the model for Galactic diffuse emission; the character c is appended to 
their names. An additional 572 (18.9%) are flagged in the catalog for less serious concerns, 
e.g., for the spectral model having a poor £t or for being close to a brighter source. Of 
the 3033 sources in the catalog, 238 (7.8%) are considered identified, based on correlated 
variability or (for 25 of the identified sources) correlated angular sizes with observations 
at other wavelengths. Of the remainder, we find likely lower-energy counterparts for 1786 
sources (59.6%). The remaining 992 sources (32.7%) are unassociated. 

The identified and associated sources in the 3FGL catalog include many Galactic and 
extragalactic source classes. The largest Galactic source class continues to be pulsars, with 
143 known 7 -ray pulsars and 24 candidates. Other Galactic source classes have continued to 
grow; fifteen globular clusters are now associated with LAT sources. Our analysis of Galactic 
source counts, informed by a model for the luminosity function, suggests that at 1 GeV ~3% 
of the Galactic diffuse emission is due to unresolved Galactic sources. Blazars remain the 
largest class of extragalactic source, with more than 1100 identified or associated with BL 
Lac or FSRQ active galaxies. Non-blazar classes of active galaxies are also found, including 
a Seyfert galaxy (Gircinus galaxy), a compact steep spectrum radio source (3G 286) and 
several radio galaxies. The populations of active galaxies in 3FGL are considered in more 
detail in the companion 3LAG catalog. 


The Fermi-LAT Gollaboration acknowledges generous ongoing support from a number 
of agencies and institutes that have supported both the development and the operation of the 
LAT as well as scientific data analysis. These include the National Aeronautics and Space 
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A. Diffuse model adjustments 

In § 13.21 we noted that the diffnse emission model has three free parameters in each 
Rol. We report their valnes in the ROIs extension of the catalog (Table [T^ and we show 
in Fignre 125] how they vary over the sky. The hrst thing to notice is that the amplitnde of 
the variations is relatively small. Overall the Galactic normalization does not vary by more 
than 20%, and in the Galactic plane (i.e., where it is the dominant component) it does not 
vary by more than 10%. The slope of the power-law correction does not exceed 0.1 (positive 
or negative) and inside the plane it does not exceed 0.05. The isotropic normalization does 
not vary by more than 40%, and ontside the Galactic plane (i.e., where it is the dominant 
component) it does not vary by more than 20%. This indicates that the diffnse model is qnite 


This preprint was prepared with the AAS LMRX macros v5.2. 


Table 15. LAT 3FGL FITS Format: ROIs Extension 


Column 

Format 

Unit 

Description 

ROI_num 

I 


ROI number (cross-reference to main table) 

RAJ2000 

E 

deg 

Right Ascension of ROI center 

DEJ2000 

E 

deg 

Declination of ROI center 

GLON 

E 

deg 

Galactic Longitude of ROI center 

GLAT 

E 

deg 

Galactic Latitude of ROI center 

Radius 

E 

deg 

ROI radius (unbinned mode) or half-side (binned mode) 

PARNAMi'"' 

E 


Value of diffuse model parameter i 

Unc.PARNAMi 

E 


la error on PARNAMi 


^Two columns (value and error) for each diffuse model parameter. The parameter name 
is given by the PARNAMi keyword in the extension header. 
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Fig. 25.— Diffuse model parameters in each Rol. The horizontal error bar is the Rol radius. 
The vertical error bar is the statistical error from the £t. The vertical scale is the same in 
the left and right plots for a given parameter. Top: Normalization of the Galactic diffuse 
component (at 500 MeV). Center: Spectral index of the power-law correction to the Galactic 
diffuse component (positive means a harder model). Bottom: Normalization of the isotropic 
component. Left: All values as a function of Galactic latitude. Right: Galactic plane only 
as a function of Galactic longitude. 
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accurate. Nevertheless, the statistical precision of the data is so good that the deviations are 
formally very signihcant. Leaving those parameters free allows releasing some of the tension 
that exists locally between the data and the diffuse model. 

At high latitudes (left-hand plots) the isotropic component is stable (as it should) but 
the Galactic normalization shows a clear North-South effect. The model is too high in the 
South but lacks emission in the North. The middle plot also indicates that the model is 
somewhat too hard particularly in the North. The group of points where the model is too 
hard in the South is around (RA,Dec) = (-M0°,—60°). 

At low latitudes (right-hand plots) the error bars on the Galactic diffuse parameters are 
very small and the parameter values are very correlated between an Rol and its neighbor. 
This is because the distance between the centers of neighboring Rols in the plane (a few 
degrees) is much smaller than their diameter (15 to 20°). The model appears to be too low 
and somewhat too hard East of the Galactic center (around longitude -|-35°) whereas it is 
too high and too soft just West of the Garina region (around longitude —80°). Outside those 
two regions the htted Galactic model is very close to the original one. Inside the plane the 
isotropic component is a minor contributor and it tends to fluctuate a lot. 

Figure [26] is another illustration of the same effect. It shows the residuals between 
the data and the full sky model (original diffuse model -|- 3FGL sources, without any free 
parameter), integrated from 100 MeV to 100 GeV. The hgure is restricted to the Galactic 
plane because at the scale shown here (0?5 pixels) nothing comes out clearly at high latitude. 
The units are sigma units (statistical deviations). The model was originally computed with 
0?1 pixels in order to model the sources accurately. The main features are the data excess 
East of the Galactic center and the data dehcit around —80°, which correspond to the 
features in Figure 1251 


B. Description of the FITS Version of the 3FGL Catalog 

The FITS format version of the 3FGL catalo^^has four binary table extensions. The ex¬ 
tension LAT_Point_Source_Catalog Extension has all of the information about the sources, 
including the monthly light curves fTable [T6|) . 

The extension Hist_Start lists the Mission Elapsed Time (seconds since 00:00 UTG on 
2000 January 1) of the start of each bin of the monthly light curves. The hnal entry is the 
ending time of the last bin. 


^®The file is available from the Fermi Science Support Center, http://fermi.gsfc.nasa.gov/ssc 
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Fig. 26.— Residuals when setting the diffuse model normalizations to 1 and no power-law 
correction, integrated from 100 MeV to 100 GeV and expressed in sigma units over 0?5 pixels. 
Top: Positive Galactic longitudes from the anticenter to Gygnus. Genter: Galactic ridge. 
Bottom: Negative Galactic longitudes from Garina to the anticenter. 
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The extension GTI is a standard Good-Time Interval listing the precise time intervals 
(start and stop in MET) inclnded in the data analysis. The nnmber of intervals is fairly 
large because on most orbits (~95 min) Fermi passes through the South Atlantic Anomaly 
(SAA), and science data taking is stopped during these times. In addition, data taking is 
briefly interrupted on each non-SAA-crossing orbit, as Fermi crosses the ascending node. 
Filtering of time intervals with large rocking angles, other data gaps, or operation in non¬ 
standard configurations introduces some more entries. The GTI is provided for reference 
and would be useful, e.g., for reconstructing the precise data set that was used for the 2FGL 
analysis. 

The extension ExtendedSources (format unchanged since 2FGL) contains information 
about the 25 spatially extended sources that are modeled in the 3FGL catalog, including 
locations and shapes. 
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Table 16. LAT 3FGL FITS Format: LAT_Point_Source_Catalog Extension 


Column 

Format 

Unit 

Description 

Source_Name 

18A 


Official source name 3FGL JHHMM.m+DDMM 

RAJ2000 

E 

deg 

Right Ascension 

DEJ2000 

E 

deg 

Declination 

GLON 

E 

deg 

Galactic Longitude 

GLAT 

E 

deg 

Galactic Latitude 

Conf_68_SemiMajor 

E 

deg 

Long radius of error ellipse at 68% confidence 

Conf_68_SemiMinor 

E 

deg 

Short radius of error ellipse at 68% confidence 

Conf_68_PosAng 

E 

deg 

Position angle of the 68% long axis from celestial North, 
positive toward increasing RA (eastward) 

Conf_95_SemiMajor 

E 

deg 

Long radius of error ellipse at 95% confidence 

Conf_95_SemiMinor 

E 

deg 

Short radius of error ellipse at 95% confidence 

Conf_95_PosAng 

E 

deg 

Position angle of the 95% long axis from celestial North, 
positive toward increasing RA (eastward) 

ROI_num 

I 


ROI number (cross-reference to ROIs extension) 

Signif_Avg 

E 


Source significance in cr units (derived from Test Statistic) 
over the 100 MeV to 300 GeV band 

Pivot_Energy 

E 

MeV 

Energy at which error on differential flux is minimal 

Flux_Density 

E 

cm~^ MeV“^ 

Differential flux at Pivot_Energy 

Unc_Flux_Density 

E 

cm~^ MeV“^ 

Icr error on differential flux at Pivot_Energy 

SpectraLIndex 

E 


Best fit photon number power-law index: for LogParabola spectra, 
index at Pivot_Energy; for PL(Super)ExpGutoff spectra, low-energy index 

Unc_SpectraLIndex 

E 


1(7 error on SpectraLIndex 

FluxlOOO 

E 

cm s ^ 

Integral photon flux from 1 to 100 GeV 

Unc_FluxlOOO 

E 

— 9 —1 

cm s 

1(7 error on integral photon flux from 1 to 100 GeV 

Energy _FluxlOO 

E 

erg cm~^ 

Energy flux from 100 MeV to 100 GeV obtained by spectral fitting 

Unc_Energy_FluxlOO 

E 

erg cm~^ 

1(7 error on energy flux from 100 MeV to 100 GeV 

Signif.Curve 

E 


Significance (in a units) of the fit improvement between power-law 

and either LogParabola (for ordinary sources) or PLExpGutoff (for pulsars) 

A value greater than 4 indicates significant curvature 

SpectrumType 

18A 


Spectral type (PowerLaw, LogParabola, PLExpGutoff, PLSuperExpGutoff) 

beta 

E 


Gurvature parameter (0 of Eq. ^ for LogParabola; NULL for other spectral types 

Unc.beta 

E 


1(7 error on /3 for LogParabola; NULL for other spectral types 

Gutoff 

E 

MeV 

Gutoff energy {Ec of Eq. [2]) for PL(Super)ExpGutoff; NULL for other spectral types 

Unc_Cutoff 

E 

MeV 

1(7 error on cutoff energy for PL(Super)ExpGutoff; NULL for other spectral types 

Exp_Index 

E 


Exponential index (6 of Eq. [2j| for PLSuperExpGutoff; NULL for other spectral types 

Unc_Exp_Index 

E 


1(7 error on exponential index for PLSuperExpGutoff; NULL for other spectral types 

PowerLaw_Index 

E 


Best fit power-law index; equal to SpectraLIndex if SpectrumType is PowerLaw 

FluxSO.lOO 

E 

cm ^ s ^ 

Integral photon flux from 30 to 100 MeV (not filled) 

Unc_Flux30.100 

2E 

— 9 — 1 

cm s 

1(7 lower and upper error on integral photon flux from 30 to 100 MeV (not filled) 

nuFnuSO.lOO 

E 

erg cm“^ 

Spectral energy distribution between 30 and 100 MeV (not filled) 

Sqrt.TS30.100 

E 


Square root of the Test Statistic between 30 and 100 MeV (not filled) 

FluxlOO.300 

E 

— 9 —1 

cm s 

Integral photon flux from 100 to 300 MeV 

Unc.Fluxl00.300 

2E 

— 9 — 1 

cm s 

1(7 lower and upper error on integral photon flux from 100 to 300 MeV^ 

nuFnulOO.300 

E 

erg cm“^ 

Spectral energy distribution between 100 and 300 MeV 

Sqrt.TS100.300 

E 


Square root of the Test Statistic between 100 and 300 MeV 

FluxSOO.lOOO 

E 

— 9 —1 

cm s 

Integral photon flux from 300 MeV to 1 GeV 
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Table 16—Continued 


Column 

Format 

Unit 

Description 

Unc_Flux300.1000 

2E 

— 9 —1 

cm s 

Icr lower and upper error on integral photon flux from 300 MeV to 1 GeV^ 

nuFnuSOO.lOOO 

E 

erg cm~^ s~^ 

Spectral energy distribution between 300 MeV and 1 GeV 

Sqrt.TS300.1000 

E 


Square root of the Test Statistic between 300 MeV and 1 GeV 

FluxlOOO.3000 

E 

— 9 —1 

cm s 

Integral photon flux from 1 to 3 GeV 

Unc.Fluxl000.3000 

2E 

— 9 —1 

cm s 

Icr lower and upper error on integral photon flux from 1 to 3 GeV^ 

nuFnulOOO.3000 

E 

erg cm~^ s~^ 

Spectral energy distribution between 1 and 3 GeV 

Sqrt.TS1000.3000 

E 


Square root of the Test Statistic between 1 and 3 GeV 

FluxSOOO.lOOOO 

E 

cm ^ s ^ 

Integral photon flux from 3 to 10 GeV 

Unc.Flux3000.10000 

2E 

— 9 —1 

cm s 

Icr lower and upper error on integral photon flux from 3 to 10 GeV^ 

nuFnu3000.10000 

E 

erg cm“^ s~^ 

Spectral energy distribution between 3 and 10 GeV 

Sqrt.TS3000.10000 

E 


Square root of the Test Statistic between 3 and 10 GeV 

FluxlOOOO.lOOOOO 

E 

cm ^ s ^ 

Integral photon flux from 10 to 100 GeV 

Unc.FluxlOOOO.100000 

2E 

cm“^ 

Icr lower and upper error on integral photon flux from 10 to 100 GeV^ 

nuFnulOOOO.lOOOOO 

E 

erg cm“^ s~^ 

Spectral energy distribution between 10 and 100 GeV 

Sqrt.TSlOOOO.lOOOOO 

E 


Square root of the Test Statistic between 10 and 100 GeV 

Variability .Index 

E 


Sum of 2x log(Likelihood) difference between the flux fitted in each time 
interval and the average flux over the full catalog interval; a value greater 
than 72.44 over 48 intervals indicates <1% chance of being a steady source 

Signif.Peak 

E 


Source significance in peak interval in cr units 

Flux.Peak 

E 

cm s ^ 

Peak integral photon flux from 100 MeV to 100 GeV 

Unc.FluxT’eak 

E 

cm~^ 

Icr error on peak integral photon flux 

TimeJ^eak 

D 

s (MET) 

Time of center of interval in which peak flux was measured 

Peak_Interval 

E 

s 

Length of interval in which peak flux was measured 

Flux_History 

48E 

— 9 —1 

cm s 

Integral photon flux from 100 MeV to 100 GeV in each interval (best fit from 
likelihood analysis with spectral shape fixed to that obtained over full interval) 

Unc_FluxJIistory 

2 X 48E 

cm~^ s“l 

Icr lower and upper error on integral photon flux in each interval 
added in quadrature with 2% systematic component 

Extended_Source_Name 

18A 


Cross-reference to the ExtendedSources extension for extended sources, if any 

0FGL_Name 

18A 


Name of corresponding OFGL source, if any 

lFGL.Name 

18A 


Name of corresponding IFGL source, if any 

2FGL_Name 

18A 


Name of corresponding 2FGL source, if any 

lFHL_Name 

18A 


Name of corresponding IFHL source, if any 

ASSOC.GAMl 

18A 


Name of likely corresponding lAGL source 

ASSOC.GAM2 

18A 


Name of likely corresponding 3EG source 

ASSOC.GAM3 

18A 


Name of likely corresponding EGR source 

TEVCATJ^’LAG 

A 


P if positional association with non-extended source in TeVCat 

E if associated with a more extended source in TeVCat, N if no TeV association 

ASSOC.TEV 

24A 


Name of likely corresponding TeV source from TeVCat 

CLASSl 

5A 


Class designation for associated source; see Tabled 

ASSOCl 

26A 


Name of identified or likely associated source 

ASSOC2 

26A 


Alternate name of identified or likely associated source 

Flags 

I 


Source flags (binary coding as in Table 


^Separate Icr errors are computed from the likelihood profile toward lower and larger fluxes. The lower error is set equal to NULL 
and the upper error is derived from a Bayesian upper limit if the Icr interval contains 0 {TS < 1, see § fT5ll . 

'^Each condition is indicated by one bit among the 16 bits forming Flags. The bit is raised (set to 1) in the dubious case, so that 
sources without any warning sign have Flags = 0. 









